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INTRODUCTION 

The baseline communications link for the Space Shuttle involves two way transmission of 
digital information via a relay satellite system. These links as presently planned appear 
to be marginal. The purpose of this study was to evaluate the feasibility of using data 
compression to improve link efficiency as an alternative to increased transmitter power, 
reducing receiver noise figures, increasing antenna gain through more stringent Orbiter 
attitude constrains, etc. 

Delco Electronics has evolved a method of encoding digital data which permits low band- 
width encoding as well as a unique system of adaptive run length encoding. The purpose 
of this study was to evaluate the effectivity of these techniques for the air-to-ground link 
and for the bandwidth-limited ground-to-ground data link used for the Orbiter downlink data. 

2.0 SYNOPSIS OF TASKS 

2.1 TASK 3. 2.1 REDUNDANCY REMOVAL, OR BITER-TO-GROUND LINK 

2.1.1 Purpose and Scope 

The purpose of this task is to establish the feasibility of improving the overall link perfor- 
mance by using data compression to reduce the transmitted bit rate. Conventional methods 
for data compression require that a substantial number of overhead bits be added to the 
transmitted data; for time division multiplexed telemetry data, the overhead management 
is such that negative throughput gains are frequently produced. The Delco system provides 
an overhead management technique which assures that all overhead bits added to the trans- 
mitted data result in a net decrease in the quantity of bits transmitted. 

The baseband Delco data, however, is not conventional in character; the baseband signal 
is a two-voltage, digital waveform composed of transitions between voltage levels in incre- 
ments of 1.5 to 4.5 bit times in one -half bit intervals. Within this family of pulses, there 
exist several which are not used to encode serial digital source data. These "unique 
characters" are then used to signal both the existence of and the quantity of redundancy. 

To evaluate the effectivity of the data compression system it is necessary to determine 
the performance of the Delco data in noisy channels. 
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The primary thrust of this task is, therefore: 

• Compare the bit error rate versus signal-to-noise ratio of the Delco 
system with that of a well known modulation technique to determine 
the compression required 

• Determine the throu^put obtainable after compression as a function 
of redundancy content 

• Estimate the redundancy content of the downlink telemetry data 

• Compare the compression ratio required with the throughput obtain- 
able to determine feasibility 

• Generate a preliminaiy mechanization incorporating the compression 
system in the air-to-ground link. 

2.1,2 Significant Results Obtained 
2. 1.2.1 Bit Error Rate Comparison 

In order to compare the Delco system with (optimum) systems emplo 3 dng conventional 
modulation techniques, optimum Delco Level I systems have been identified and an analysis 
of their performance carried out. This subtask has been accomplished and is reported 
in Appendix B. The study has been restricted to baseband systems, but comparisons among 
baseband systems generally carry over to the corresponding modulation systems. (See Ref- 
erence 1 of Appendix B. ) Optimum Level I systems have been identified and studied for 
both the usual Level I application in which there is a bandwidth constraint and for appli- 
cations that have no significant bandwidth restrictions. These systems are not truly 
optimum in the sense that their performance can be improved by utilizing sequential de- 
coding techniques due to the correlation in the Level I bit streams. However, sequential 
detection of Level I has not been attempted in the laboratory, and its analysis was considered 
outside the scope of the study. Sequential detection of Level I will improve the error per- 
formance at the cost of some increase in hardware complexity. 

The basic results of the study are shown in a simplified form in Table 1. For more 
detailed results and discussion, see Figure 23 and Section 4,4 of Appendix B. 


MODULATION TECHNIQUE 

SNR TO ACHIEVE A lO"^ BER 

Conventional 

8.4 dB 

Conventional with Coding Gain 

3.8 dB 

Level I 

11.7 dB 

Level I (Bandwidth Restricted) 

13.5 dB 
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The signal-to~noise ratio (SNE) is defined to be the ratio of the received signal power to 
the noise power in a bandwidth equal to the information rate. The entry labeled "Conven- 
tional with Coding Gain" is extracted from the article by Batson and Moorehead ("A Digital 
Communications System for Manned ^aceflight Application," B. H, Batson and R.W. Moore- 
head) in which a 4. 6 dB reduction in required signal-to-noise ratio is shown with the use 
of convolutional encoding and Viterbi decoding. 

Table 2 shows the reduction in signal-to-noise ratio requirements that would be realized 
through the use of Level II and Level III as a function of the compression ratio obtained. 


Compression Ratio 

1:1 

2.1:1 

6.2:1 

12.3:1 

Required SNR (before 

11.7 dB 

8.4 dB 

3.8 dB 

0.8 dB 

compression to achieve 





11. 7 dB SNR after compression) 
1 





^4 

Table 2. SNR to Achieve 10 BER Versus Compression Ratio 


This table has been constructed from the equation (SNR reduction) = 10 log (compression 
ratio) and, of course, is valid regardless of the means used to obtain the compression 
and thus reduce the data rate. The entries were chosen to show that a compression of 
2.1:1 suffices to bring I^vel I up to the performance of the conventional system; a com- 
pression of 6.2:1 achieves the error rate of the conventional system with convolutional- 
Viterbi; and a compression of 12.3:1 reduces the SNR requirements to 3 dB below those 
of coded conventional. 

2 . 1 . 2 . 2 Throughput Improvement 

The compression system considered is one in which the sampled signal value is sent during 
selected master PCM frames or when the signal values differ from a previous sample by 
a predetermined amount. When a signal or string of signals are not sudfficiently different 
from previously transmitted values, uniquely coded waveforms are substituted for these 
signals. These waveforms occupy significantly less space in the communications channel 
and permit a throughput improvement. Buffering of the information to be transmitted 
provides for the generation of a constant bit rate PCM serial output, which carries aU in- 
formation contained in the original PCM data, but at a much lower bit rate. The reduc- 
tion in bit rate achievable is related to the redundancy contained in the source data and to 
the overhead required to identify the data channels which were not sent in original signal form. 
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This system generates no overhead during nonredundant periods and from 4 to 23 bits for 
redundant strings of from 1 to 256 data channels* The compression produced varies from 
1:1 for zero redundancy to over 90:1 for string lengths of approximately 250. The system 
adaptively selects the preferred method for encoding the output PCM data — for nonredundant 
channels, the data are encoded in Level I; for redundant strings of up to eight bits, the data 
are encoded by Level II; for strings of nine bits or more, the data are encoded by Level HI 
coding. The mode of operation selected is easily detected by the receiving system due to 
the difference in the use of the unique characters used to encode the redundancy. 

The differences between the encoding levels can be seen from examination of Figures 1 
through 4. The rules for encoding data in Level I are portrayed in Figure 1. In this 
mode of operation, the unique character is used to encode three consecutive encoded 
zeros; this eliminates the possibility of generating large, low frequency components in 
the baseband Level I signal. When the imique character is used for this purpose, it is 
injected as near the preceeding transition as is permitted by the encoding rules. It is 
never delayed from this position when used to encode data by Level L 

Figure 2 illustrates how this character is delayed to produce Levels II and HI; the two 
following figures are examples of that encoding. The Level III mode of operation is indi- 
cated by two sequential, unique characters, each with a one bit injection delay. The run 
length counter is a binary count of the run length, transmitted most significant bit first. 

The counter value is always followed by a bit controlled to a logic ONE state and a delayed 
unique character. This arrangement clearly indicates the presence of Level III encoding, 
and the coimter length is indicated by the terminator coding. 

Figure 5 illustrates the relationship between string length and overhead and the achieve- 
able compression, 

2 . 1 . 2 . 3 Downlink Redundancy Content 

In order to arrive at an achievable througl^ut improvement for Orbiter data, a sample 
PCM format was constructed using selected Titan data signals in expanded quantities, 

A 360-word format was synthesized to establish the serial stringing produced by the re- 
dundancy contained in the sampled signal values. This format is contained in Appendix A. 
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Reference 

Transition 


\ I 


One Redundant Channel; 4 Bits Overhead; 4 Bits Saved 



Two Redundant Channels; 5 Bits Overhead; 11 Bits Saved 



Three Redundant Channels; 9 Bits Overhead; 15 Bits Saved 



T "1 I Four Redundant Channels; 10 Bits Overhead; 22 Bits Saved 



"i r ”• T Five Redundant Channels; 11 Bits Overhead; 29 Bits Saved 



NOTE 1. Bits Saved are Based on 8- Bit PCM Data 
NOTE 2 . Reference Transition is First Transition 
in Word Preceeding Redundant Channels; 
Remainder Sent Following Redundant Word Coding 


{ I Eight Redundant Channels; 17 Bits Overhead; 

47 Bits Saved 


Figure 3. Level II Encoding 


Source Encoding of Data That Includes Short and Long Redundant Data Strings 

LEVEL Ml 


Reference 

Transition 


LEVEL 1 1 


_l L 







One Redundant Block 


1 Delay + 1 Delay = Run Length Run Length 
Start of Run-Length Code of Terminator. 

Code 110 1 

Mode of Operation is Adaptively Selected, For Run Lengths of 7 or 
Less Level 1 1 is Used. For Runs of Over 8, Level III is Used. 


• Run Length Coded MSB First; Number of Bits 
Used is Variable and Dependent on Redundancy 
Present. MSB is Always a One and Immediately 
Follows Start Code. 


Figure 4. Level III Encoding 
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Figure 5. String Length/Overhead Relationships 
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These signals were examined during the boost and coast phases, the signal activity during 
the coast period was selected as being most representative of Orbiter on-orbit data condi- 
tions where the weak link conditions exist. During this period the signals are highly redun- 
dant, and fewer than two data channels require transmission out of each frame of data. A 
throughput advantage of over 100:1 is suggested; however, periodic signal updates, synch- 
ronization codes, and a reasonable safety factor reduce this factor significantly. A more 
reasonable design goal is 10 to 20:1, which permits a greater variation to exist due to signal 
anomalies without overloading the compression system. The redundancy present in the 
boost phase, where maximum signal activity exists, permits compression of approximately 
10:1. The string length analysis for those flight phases are also contained in Appendix A. 

2. 1.2.4 Feasibility 

The results of Section 2. 1.2.1 above show that if Levels I, n, and III can develop a com- 
pression of 10:1 or better, the signal-to-noise requirements will be reduced substantially 
compared with a conventional system using convolutional encoding — Viterbi decoding. 

Section 2. 1.2* 3 conservatively estimates that compression ratios of 10 to 20:1 are achieve- 
able with the downlink data. It is therefore concluded that not only is the application of 
Delco's source encoding techniques to the Orbiter-to-ground communications links feasible, 
but that the use of Levels I, II, and ni provide a sufficient reduction in transmitted bit rate 
so as to negate the requirement for channel encoding. 

2 . 1 . 2 . 5 Me chanization 

A generalized data compression system is shown in Figures 6 and 7. This system 
accumulates data from the existing PCM Master Unit (PCMMU). Each data channel is 
compared with its previous value, and channels which contain new information are tagged 
for transmission. The buffered data are searched for these tags in an order identical to 
the original PCM format. All channels tagged are included in the output data stream. It is 
apparent that a variable word/frame rate exists for the output data. The variable word/ 
frame rate is restored by buffering in the ground system for data playback in coiventional 
data systems, if required. 
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•1 MEMORY LOCATION FOR EACH 
TELEMETERED VVORD (8 BIT) 


•11 BITS SEMICONDUCTOR RAM PLUS 


4 BITS SEMICONDUCTOR PROM FOR 
EACH LOCATION 


DATA WORD 


TRANSMIT 

CONTROL 


FILTER 

COUNT 



APERTURE SIZE 


<£> 

Figure 6. Airborne Compressor Functional Diagram 
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Figure 7, Ground Based Decompression System Functional Diagram 
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2.1.3 Conclusions 

The telemetry portion of the Orbiter air/ground link can be compressed sufficiently to pro- 
vide a significant increase in channel performance. The voice portion of the data link must 
be similarly compressed in order to achieve the desired performance increase. It is 
understood that NASA has separately contracted for a study concerned with voice compression. 

2.1.4 Recommended Action 

It is recommended that followon studies be initiated to permit an evaluation of the per- 
formance of Level I when combined with spread spectrum. It is anticipated that the longer 
transition intervals guaranteed by the Level I data will produce meaningful processing gains 
and enhance synchronization techniques. 


2.2 TASK 3.2.2 REDUNDANCY REMOVAL FOR ORBITER GROUND-TO-GROUND LINK 

2.2.1 Purpose and Scope 

The purpose of this primary task is to determine the potential for using data compression 
techniques to provide improved ground link data throughput and to consider methods for 
improving the integrity of the transmitted data signal through the use of Delco three-phase 
encoding. Unlike the air-to-ground data link, these channels are bandwidth limited, and 
data compression is desirable to decrease the cost of multiple, wideband data links, which 
would be necessary to convey the uncompressed source data. In order to accomplish the 
intent of this task, the following items were accomplished; 

• Bit error rate studies to compare Delco encoding with conventional methods 

• Estimate the redundancy content of Orbiter PCM data 

• Determine achievable throughput using data compression 

• Evaluate the compatibility of Delco encoding techniques with channel 
encoding and decoding 

• Generate a preliminary mechanization that illustrates the method of incor- 
poration of selected techniques. 

2.2.2 Significant Results Obtained 
2. 2. 2.1 Bit Error Rate Studies 
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2. 2. 2. 1. 1 Level I in Bandwidth Limited Applications 

Section 2 of Appendix B is devoted to the definition and analysis of an optimum Level I 
system for bandwidth limited applications. The optimum system identified is very 
similar to the Delco laboratory system. The two systems differ, at most, by nonessential 
differences in the transmitting and receiving filters. This is an important fact. The 
optimum Level I system transmits data at the Nyquist rate; that is, at a rate which is twice 
the bandwidth of the channel. The laboratory system also achieves transmission at the 
Nyquist rate utilizing simple, low cost hardware. Most systems operating over band- 
width-limited channels operate at only one-half to two-thirds the Nyquist rate in order to 
reduce complexity and simplify the hardware requirements. Delco' s Level I achieves full 
Nyquist rate transmission without significant reduction in error performance. 

Table 1 illustrated the basic results. Level I in bandwidth-restricted applications has only 
1.8 dB less margin over noise than fully expanded Level I. The results of Table 2 may, 
therefore, be adapted to bandwidth- restricted Level 1 by adding 1.8 dB to each of the signal- 
to-noise ratio entries there. 

2. 2. 2. 1.2 Delco Electronics 3-Phase Versus Conventional 4 -Phase Differential 
Phase Modulation 

An analytical expression was developed for phase error probability as a function of signal- 
to- noise ratio for L-phase differential phase modulation (DPM) and evaluated numerically 
for biphase, 3-phase, and 4-phase modulation at signal-to-noise ratios ranging from 0 to 
15 dB. Subsequently, the effectiveness of Delco Electronics 3-phase versus conventional 
4-phase DPM was determined, (Refer to Appendix C.) 

The results indicate that, for example, at a SNR of 10 dB, the probability of a phase error 
in 3-phase modulation is down by a factor of ten from that of 4-phase, As SNR increases, 
this factor increases, and the advantages of 3-phase modulation become more prominent. 
Translating phase error rate into bit error rate, our results for conventional 4-phase 
modulation show that a single phase error produces on the average 1 bit error; for Delco ’s 
3-phase modulation, a phase error produces on the average approximately 1-1/2 bit errors. 
Comparing the two techniques at 10 dB, the 3-phase modulation bit-error rate is down by 
a factor of about seven from that of 4-phase modulation. 
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2. 2. 2. 3 Compatibility of Convolutional and Level I Encoding 

A study was conducted to determine if a convoluted code could be encoded in Level I 
without impairing the error correcting effectiveness of Viterbi decoding. Rate one-half 
convolutional encoding of random source data was used for the investigation. 

Two digital computer programs were employed in error-correcting simulations. (Refer 
to Appendix D, Section 5.) The first permits rate one-half convolutional encoding, error 
injection, and Viterbi decoding. The second incorporates rate one-half convolutional en- 
coding followed by Level I encoding, error injection. Level I decoding (using the Viterbi 
decoding algorithm), and Viterbi decoding. 

Simulations were made with a source data sequence of 4,500 random bits. Error rates 
were established which were high enough to result in errors in the decoded sequence. 

For random source data. Level I encoding leads effectively to a bandwidth expansion by 
a factor of two. Taking this into account throu^ the expression for the probability of a 
channel error in terms of bandwidth, the results showed that the error correcting ability 
of Viterbi decoding is not decreased at hi^ error rates. Below an error rate of about 
30%, the error correcting ability of Viterbi decoding is, however, degraded through the 
introduction of Level I. 
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2 . 2 , 2 . 2 Redundancy E s timate 

2 . 2. 2 . 2. 1 Priority Transmission - Uncompressed Air/ Ground Link 

In the event that compression is not used for the downlink, priority transmission can 
provide for the real time transmission of vital data followed by nonreal time transmission 
of less critical and/or redundant information. 

Critical data consists of out-of-limits signals, changes in the state of discrete signals, and 
samples of signals selected for consumables management and other monitoring tasks. The 
percentage of signals of real time interest is very small - it has been estimated that a 
maximum of 5% of the total monitored signals will display out-of-limits conditions for 
Orbiter malfunctions, and these will persist only until redundant Orbiter systems are 
placed on-line to correct the abnormal condition. This percentage will result in approx- 
imately 800 words per second which are candidates for transmission; further processing 
to remove the redundancy in these signals will reduce the candidates significantly such that 
they can easily be accommodated by a conventional 4, 8 kb/s telephone transmission link. 
Appendix A contains the results of the redundancy reduction achieved for telemetered data 
signals, and it is reasonable to assume that out-of-limits channels will exhibit the same 
average change history as in-limits signals. Nonreal time transmission of the accumu- 
lated flight data is accomplished using a ±1, ”one" bit change criteria, which according to 
the data in Appendix A, permits transmission of the entire data set well before the next 
orbital pass. 

2. 2. 2. 2. 2 Priority Transmission - Compressed Air/Ground Link 

When compression is used in the downlink system, the degree of redundancy eliminated will 
necessarily be less than that which can safely be used in the ground link. Data eliminated 
from the downlink is generally not recoverable for post-test usage, which is not true of 
the ground link system. All data received from the downlink can be recorded for later 
use, which permits greater latitude in the ground compression system. This permits using 
redundancy detection with wider apertures and heavier filtering than would be prudent for 
the downlink case. As in the uncompressed downlink case, limit detectors and discrete 
processors provide for real time transmission of critical data and nonreal time", which are 
compatible with standard telephone data transmission techniques. 
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2.2.3 Mechanization 

The processing subsystem selected is one which has sufficient speed to assure evaluation 
of the full 16k words/second contained in the uncompressed downlink data. These pro- 
cessors are of the type produced by Delco for several military programs including the 
TDMG and SSTC equipments, and are identical to those proposed for the processor option 
for the Space Shuttle PCMMU. (Reference Delco proposal P74-2-6-1. ) The functional 
mechanization for the ground-based compressor system is shown in Figure 8. 



• CHANGED DISCRETES 

• SELECTED DATA 


Figure 8, Ground Based Data Compressor Basic Mechanization 
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2.2.4 CONCLUSIONS 

The ability to store real time data at the receiving terminals for delayed transmission 
significantly reduces the real time throughput requirement if sufficient signal processing 
is provided to extract signals of immediate interest. High speed telemetry data preproc- 
essors exist for this purpose; those proposed for the PCMMU preprocessors are examples 
of proven techniques. Use of these preprocessors for selection of out-of-limits channels, 
changes in discrete signals, and other signals of interest provides for a very low real time 
throughput requirement. For this application, the throughput rate necessary to support 
the program can be accomplished without the use of Level I or Delco 3-phase DPSK, These 
encoding systems have been shown to possess very desirable characteristics for telephone 
data transmission; however, completely adequate transmission rates can be provided for 
compressed data using conventional 4-phase DPSK and commercially available data MODEMS. 

2.2.5 RECOMMENDATIONS 

In view of the vital role played by the telemetry data preprocessors in the ground data link, 
it is recommended that a study be undertaken to prepare a detailed specification for the 
ground-based telemetry preprocessor unit. 
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APPENDIX A 

SAMPLE PCM FORMAT 
USING SELECTED TITAN IRC DATA SIGNALS 

1*0 INTRODUCTION 

The results of the analysis performed to determine the expected redundancy content of 
the Orbiter link are summarized herein. Actual Titan IIIC flight data was subjected to 
the data compression technique of determining the percent of data that exceeded programm- 
able apertures. As a result of this analysis, additional methods and mechanization modif- 
ications were identified and investigated, the results of which are also included herein. 

2.0 APPROACH 

Actual Titan IIIC telemetry data, selected from the 30 May 1974 launch, was used for 
the redundancy assessment tests via a Sigma 7 computer simulation. The basis for 
selection of the 33 vehicle and guidance analog measurements, provided in Table 1, was 
to provide; 

• Similarity to envisioned type of orbiter data 

• A mixture of sampling rates 

• A mixture of quiescent and highly active signals 

• Signal signatures that are dependent upon environmental conditions. 

The types of data utilized consisted of power supplies, discretes, steering/actuator func- 
tions, temperatures, pressures, accelerations, vehicle rate, and IGS bus voltages and 
currents. All measurement data types were subjected to the data compression techniques 
over 30 to 50 seconds of the following four desired environmental periods: (1) quiescent 
prelaunch (Ig), (2) maximum dynamic pressure, (3) maximum thrust (approximately 4. 5g), 
and (4) coast. 

The Sigma 7 computer program simulated the redundancy reduction technique of only trans- 
mitting data values when a programmed allowable aperture was exceeded. Upon exceeding 
the allowable aperture, an updated tolerance, with the transmitted data value equal to the 
new nominal, was derived and used for subsequent transmission tests. Various reasonable 
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SIGNAL TYPE 

SAMPLE RATE 
PER SECOND 

Roll, Pitch, Yaw Rate Gyro Outputs 

400 

Pitch, Yaw Lateral Accelerometer Outputs 

400 

IGS Voltage and Current 

20 

SA 1, 2, and 3 Thrust Chamber Pressures 

400 

IL 1, 2, and 3 Current Ladder Outputs 

100 

ELI and 2 Voltage Ladder Outputs 

100 

BL 2 Sequence System Discretes 

100 

BL 6 and 7 ACS Nozzle Discretes 

100 

Guidance Truss (MGC) Temperature 

100 

Guidance Truss (IMU) Temperature 

20 

Four ACS Nozzle Pressures 

200 

X, Y, Z Payload Accelerometers 

400 

IMU Internal Temperature 

20 

IMU “15 Vdc Unregulated 

20 

IMU +15 Vdc Regulated 

20 

IMU 2 8 Vac, 0 Phase 

20 

MGC -6 Vdc 

20 

TPS Bus Current 

800 


Table 1. Signals Tested for Redundancy 


magnitudes were evaluated to aid in determining the minimum bit rate required to support 
the transmission of orbiter data. 

3.0 DISCUSSION 

The results of subjecting the actual Titan IIIC flight data to the data compression technique 
of determining percent of data that exceeded programmable apertures is provided in 
Tables 2 throu^ 7. The Sigma 7 computer output was grouped into the categories listed 
on each of the tables. The group average percent per flight period is summarized for each 
aperture tested. 


Note: Tables 2 through 13 and Figures 1 and 2 appear at the end of this appendix. 
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The differences in signal signatures is evident in comparing the summarized group 
averages of active signals, provided in Tables 2 and 3, with those of very quiescent 
signals, provided in Tables 6 and 7. Hardcopy chart recordings of each of the signals 
were made to aid in the Sigma 7 program checkout, and an example of some of the active 
signals during the maximum dynamic pressure period is provided in Figure 1* 

Examination of the results initially tabulated illustrated that several signals exceeded the 
± 1 bit aperture for a significant percentage of the time, while the actual signal only ex- 
ceeded ± 1 bit less than 2% of the time. This was due to the fact that the initially stored 
signal value was at a peak value rather than the true nominal value. This problem was 
alleviated by programming additional initialization criteria, where the nominal for the 
first aperture was derived from the average of the samples contained in the first record. 


The importance of electrically terminating defunct signals when incorporating a data com- 
pression scheme such as that tested was clearly illustrated in several of the 33 signals; 
these signals are footnoted on the tables by an asterisk. These conditions occurred during 
the parking orbit coast phase with detached Stage H signals, including a vehicle roll rate 
gyro and a thrust chamber pressure. The variations exceeded the allowable ± 3 counts 
99. 8% of the time for one signal and 75% of the time for the other signal. Obviously, 
conditions like these would severely defeat data compression objectives. 


The total percentage of all 33 signals exceeding the tested apertures per flight phase was 
as follows. 



± 1 Bit 

± 2 Bits 

± 3 Bits 

+ 5 E 

Prelaunch 

8.8 

6, 6 

5.5 

4.0 

Maximum Q 

24,1 

19.2 

16.3 

12.1 

Maximum Thrust 

14,5 

7.3 

4.5 

2.2 

Coast* 

16.1 

10.3 

10.2 

10.0 

Coast 

5.7 

0.3 

0.2 

0.1 


* Includes percentages of noisy defunct signals. 


As is apparent when reviewing the group average percentages per flight phase for each of 
the categorized six groups, primary emphasis should be placed on the most active signals. 
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Which are summarized in Tables 2 and 3. Thus, additional logic was defined and imple- 
mented in the Sigma 7 simulation to define exceeding a programmable aperture only when 
this occurred for a programmable number of consecutive samples. Tables 8 and 9 contain 
the results of the following three tests: (1) greater than 1 bit for two consecutive samples, 
(2) greater than 2 bits for two consecutive samples, and (3) greater than 1 bit for three con- 
secutive samples. Direct comparison can be made between Tables 2 and 8 and between 
Tables 3 and 9. The average percentage of the four groups of active IGS data exceeding 
the ± 1 bit aperture (Table 2) was 48% compared with 22% of the same data exceeding the 
± 1 bit for two consecutive samples (Table 8), Likewise, the reduction for the i 2 bit 
is 27% compared with 9%, The percentage of data exceeding the greater than 1 bit for 
three consecutive samples is comparable to that of the greater than 2 bit for two consec- 
utive samples. Similar reductions are apparent for the active airframe measurements, 
where 48% and 38% (Table 3) are reduced to 22% and 18% (Table 9) for the ±1 and ±2 bit 
cases, respectively. 

The significant advantages of incorporating the consecutive sample logic is further sup- 
ported by the fact that the noise outputs of the defunct signals during coast (footnoted by 
an asterisk) are suppressed. Although the ± 5 bit aperture without consecutive logic would 
provide a similar data compression ratio as the ± 2 bit with consecutive logic, important 
signal characteristics could be obscured if the signal sensitivities were optimally defined. 

Mean and standard deviations of the time between aperture excursions were also computed 
for all 33 signals. However, the usefulness of the statistical outputs was very limited. 
Significance of the mean tests were only satisfied for the very active IGS and airframe 
measurements. These statistics for the active IGS measurements are provided in Table 
10, When no aperture excursions were recognized, the mean is equal to the sample size, 
and the standard deviation is zero. Also, the mean has been rounded and the standard 
deviation truncated. 

4.0 ADDITIONAL INVESTIGATIONS 

The average string length of zeros, documented herein, was calculated to determine a 
compression ratio that would be realized from a data set equal to the 33 signals tested 
at the Titan IIIC telemetry sampling rates. 
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The format of the repeatable block of 353 samples 


■ 33 

1=1 


is provided in Table 11. The percentage of data transmitted was that resulting from ex- 
ceeding the recommended greater than 2 bits for two consecutive samples data compression 
aperture. 


A very large compression ratio would be realized during coast since the only data exceed- 
ing aperture was 23% and 6% of two 20 S/s signals, 1% of two 100 S/s signals, and 1% of 
one 200 S/s signal. With these very few aperture excursions, one sample of one signal 
would be transmitted only 49% of the time. Thus, for the worst case, 51% of the time 
the string length of zeros would be the total block of 360 words. With random distribution, 
the average string length during coast for the 360-word block would be approximately 340 
words, which would yield a compression ratio of 100:1. 

To determine the achieveable compression for more active periods, the percentage of data 
exceeding the subject aperture during the period of maximum dynamic pressure (max. Q) 
was obtained. An illustration of the transmitted data in the 360-word format is provided 
in Table 12. The calculations of the average string length during maximum Q are tabulated 
in Table 13, and the result is plotted in Figure 2. As illustrated, the average string 
length leveled out at 18. 85 after a sample size greater than 80. This average string length 
would result in an average compression ratio of 10:1. 

5.0. CONCLUSION 

Dependent upon the degree of data compression desired, either the greater than 1 bit or 
greater than 2 bits for two consecutive samples should be implemented in the transmission 
of Orbiter data. The greater than 2 bits for two consecutive samples data compression 
technique is recommended to provide the most meaningful data at a highly desirable com- 
pression ratio. Ratioing these resulting percentages of the data compression technique 
selected by the number, type, and sample rate of measurements planned for Orbiter data 
will determine the predicted data compression and required bit rate for transmission of 
Orbiter data. 
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pMASr 

SIGNAL 
PIATE; #; 



APpRTU^P / 

PERCr^JT 


T0TAL- SAMPLE? 



NAME 


+ /- 

1 

^IT 

+/- 

2 

BITS 

4/ 

«. 3 

BITS 

4 /- 

5 

9TTS 

FREL;U^;C'^ 

2o;1180; 

603 

1 4g 

/ 

24*5 

1C7 

/ 

17.7 

48 

/ 

8*0 

1 4 

/ 

P.3 

Q 

100? 

PSC 

/ 

27*9 

181 

/ 

18.0 

85 

/ 

8.5 

20 

/ 

9,^ 


IGS Volts 


125 

/ 

25*8 

53 

/ 

11.0 

8 

/ 

1.7 

8 

/ 

1.7 

CbAFT 



1.00 

/ 

??*7 

54 

/ 

12.3 

23 

/ 

5.2 

13 

/ 

2*9 


2c; 1181; 

603 . 

447 

/ 

74*1 

421 

/ 

69.8 

372 

/ 

61.7 

191. 

/ 

31.7 


IGS Current 

1 003 

952 

/ 

94*9 

907 

/ 

90»4 

763 

/ 

76. 1 

321 

/ 

32*0 

Th'-UST 

A84 

443 

/ 

9t *5 

364 

/ 

75.2 

274 

/ 

56.6 

Xl-5 

/. 

44.4 

OCAS'" 


440 

967 

/ 

6C*7 

191 

/ 

43.4 

138 

/ 

31.4 

IPC 

/ 

27.0 

phEL AUNC w 

20;322; 

603 

?30 

/ 

38*1 

1^ 

/ 

2*3 

0 

/ 

*0 

C 

/ 

.0 

t^AxI^^U^ Q 


3 003 

501 

/ 

50*0 

2 

/ 

.2 

0 

/ 

.0 

G 

/ 

• 0 

TM^UST 

MGC-6Vdc 

484 

241 

/ 

49*8 

2? 

/ 

4 . 5 

0 

/ 

• 0 

C 

/ 

*0. 

C&aST 


440 

24 

/ 

5*5 

0 

/ 

• 0 

0 

/ 

.0 

0 

/ 

.0 

PRL t, 

2c; 228; 

603 

1B7 

/ 

26*0 

llB 

/ 

19.1 

0 

/ 

.0 

c 

/ 

• 0 

'■'iAx.T-UM Q 

Temp., IMU 

1003 

437 

/ 

4 3*6 

0 

/ 

*0 

0 

/ 

.0 

0 

/ 

.0 

^Ax* Thrust 

CbAST 

4 84 

173 

/ 

35*7 

62 

/ 

12.8 

0 

/ 

•0 

c 

/ 

. 0 

Elat. 

440 

P99 

/ 

68*0 

20^ 

/ 

46.4 

132 

/ 

30*0 

0 

/ 

.0 


Group Avg. % 
for: 


Prelaunch 

40.7 

27.2 

17.4 

8,4 

Max. Q 

54.1 

27.2 

21.2 

8.5 

Max. Thrust 

50.7 

25.9 

14.6 

11.5 

Coast 

39.2 

25.5 

9.2 

7.6 


> 

I 

C5 


Table 2. Active IGS Measurements Throughout Flight 
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/ 

PERCE'^'T ?F reiAL samples 


+/- 1 

BTT 

■«-/- 

2 

2 ITS 

+ /- 

3 

BITS 

+ / ■ 5 

3 

ITS 

2 CJ /+ Q 

/ 

2'+. 5 

2001 

/ 

16.6 

1506 

/ 

12.5 

885 

/ 

7.? 

131?6 

/ 

65.4 

97 64 

/ 

48.7 

7459 

/ 

37.2 

4332 

/ 

21.6 

1?63 

/ 

13.0 

343 

/ 

3.5 

141 

/ 

1.5 

8 

/ 

. 1 

*3? 19 

/ 

37.7 

c 

/ 

.0 

C 

/ 

.0 

0 

/ 

.0 

2722 

/ 

22 . 6 

lS6l 

/ 

15.4 

1346 

/ 

11.2 

741 

/ 

6.1 

10775 

/ 

53*7 

7333 

/ 

36.5 

5192 

/ 

25.9 

2814 

/ 

14.0 

35^0 

/ 

36.6 

2220 

/ 

22.9 

1501 

/ 

15.5 

791 

/ 

8 • 2 

0 

/ 

*0 

0 

/ 

.0 

0 

/ 

.0 

r . 

/ 

.0 

33S1 

/ 

2S*0 

•'2487 

/ 

20.6 ■- 

1980 

/ 

16.4 

1274 

/ 

10.6 

US30 

/ 

73*9 

12028 

/ 

60.0 

9765 

/ 

48. 7 

6410 

/ 

32.0 

2557 

/ 

26*4 

1113 

/ 

11.5 

553 

/ 

5.7 

130 

/ 

1.3 

♦ j? 7 79 

/ 

9?. 8 

877'=* 

/ 

99.8 

8779 

/ 

99.8 

8860 

/ 

99.8 

4347 

/ 

36.0 

3570 

/ 

29.6 

3256 

/ 

27.0 

2777 

/ 

23 • 0 

18652 

/ 

93.0 

17694 

/ 

88.2 

16794 

/ 

83.7 • 

l5085 

/ 

75.2 

6289 

/ 

65*0 

4612 

/ 

47.6 

3400 

/ 

35.1 

197C 

/ 

20.4 

0 

/ 

.0 

0 

/ 

.0 

0 

/ 

• 0 


/ 

r>. 

• \ . 

3632 

/ 

3C.1 

3l7'3 

/ 

26.4 

2845 

/ 

23.6 

230A 

/ 

19.1 

17791 

/ 

88.7 

16377 

/ 

81.6 

15141 

/ 

75.5- 

12754 

/ 

63*6 

4215 

/ 

43 . 5 

2443 

/ 

25.2 - 

154b 

/ 

15.9 •• 

608 

/ 

6.3 

0 

/ 

*0 

0 

/ 

.0 

0 

/ 

.0 

0 

/ 

.0 


28.2 

21.7 

18.1 

13.2 

74.9 

63.0 

54.2 

41.3 

36.9 

22.1 

14.7 

7.3 

27.5 

20.0 

20,0 

20.1 


Powered Flight, Active Airframe Measurements 


DCLCO ELECTRONICS DIVISION • SANTA BARBARA OPERATIONS • GENERAL MOTORS CORPORATION 



R74-34 A~8 



SIGNAL 




P A S E 

RATE; #; 

NUMBER 

«»«*«*««#« 


NAME 

Samples 

1 

PIT 


iod; 301 

3015 

66 / 

2.2 


IL 1 

5015 

PCS / 

4.2 

KAx. THRUST 



2 / 

.1 

CeAST 


2E0G 

0 / 

.0 

PRELAUNCH 

loo; 302; 

3015 

75 / 

2.5 


TT O 

5015 

165 / 

3.3 

max* thrust 

IL 2 

2430 

3 / 

.1 

CbAST 


2200 

0 / 

• 0 

PFELA'JNCH 

lOo; 303; 

3015 

62 / 

2.1 

^^AXT“UH Q 


5015 

3 36 / 

2*7 

MAX* Thrust 

IL 3 

2420 

4 / 

• 2 

CbAST 


220C 

0 / 

*0 

PREL a'JNCH 

100; 310; 

3015 

134 / 

4.4 

maximum q 

EL 1 

5015 

214 / 

4.3 

max. ThouST- 


2420 

Q / 

.0 

CbAPT 


2200 

0 / 

• 0 

prelaunch 

ioo;311; 

3015 

114 / 

3. a 

maximum Q 

EL 2 

5015 

?91 / 

5*8 

max. THPUST 


2420 

2 / 

.1 

CttAST 


2200 

0 / 

.0 

Group Avg, % for; 




Prelaunch 




3.0 

Max, Q 




4.1 

Max. Thrust 




0.1 

Coast 




0.0 


APrRTUr^r. / PR^Ci^'^'T TBTAL SA^iPLE"^ ♦**#**##** 


+/ 

• 2 

HITS 

+ /- 

3 BITS 

+/- 

5 

BITS 

5^ 

/ 

1.7 

40 

/ 

1*3 

35 

/ 

l.P 

103 

/ 

2.1 

53 

/ 

1.1 

23 

/ 

.b 

0 

/ 

.0 

0 

/ 

.0 

p 

/ 


0 

/ 

.0 

0 

/ 

.0 

0 

/ 

.0 

54 

/ 

1.8 

44 

/ 

1.5 

34 

/ 

1.1 

55 

/ 

1.1 

42 

/ 

.8 

1^- 

/ 

.3 

n 

/ 

.0 

0 

/ 

.0 

n 

/ 

• tj 

0 

/ - 

.0' 

0 

/ 

.0 


/ 

• ‘ ‘ 

51 

/ 

1.7 

45 

/ 

1.5 

36 

/ 

1.2 

67 

/ 

1.3 

31 

/ 

. 6 

14 

/ 


0 

/ 

.0 

0 

/ 

.0 

0 

/ 

♦ L 

0 

/ 

.0 

0 

/ 

.0 


/ 

.D 


/ 

3.8 

103 

/ 

3.4 

25 

/ 

<?• 

# 

370 

/ 

3.4 

135 

/ 

2.7 

17 

/ 

• 3 

0 

/ 

.0 

0 

/ 

.0 

C 

/ 

.0 

0 

/ 

.0 

.0 

/ 

.0 

0 

/ 

.0 

30 

/ 

1.0 

26 

/ 

.9 

19 

/ 

. 6 

191 

/ 

3.8 

137 

/ 

2.7 

55 

/ 

1*1 


/ 

• 1 

0 

/ 

.0 

0 

/ 

.0 

0 

/ 

.0 ■ 

0 

/ 

.0 

V ' 

/ 

.0 


2.0 

1.7 

1.0 

2.3 

1.6 

0.5 

0.0 

0,0 

0.0 

0.0 

0.0 

0.0 


Table 4, Powered Flight, Steering Ladder Outputs 
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phase 


SIGNAL 

RATE; 

NAME 


number 

Samples 


APCRTtjf^E / PERCENT 
2 HITS 


?p T0TAL SAMPLE' 
+/- 3 PITS 


»***•»****■* 
>/- 5 PFTS 


prelauncw see; 1179; 

M ' V T ^ ^ Q 

HAftt thpust "T1*S Bus Current 
CL ' A ST 


PREL *’-'^ C '-» 

M^yjMlJM Q 

max * T *- J ' US ' 

COAST 

ppel^^^^ch 
HaxIHU '* Q ' 
max * Thpijs 
coast 


aoc;2350; 

P/L Accel. 
Long. 

aoc;2351; 

P/L Accel. 

Lat. 


pr£ launch 
MAxI""UM Q 
max. TIPUST 

CuaST 


Aoci 2352 ; 

P/L Accel. , 
Vert. 


prela'Jncm AocJSOlO; 

TMf:!isTT/C Press SA3 
coast 

pp£;LiUHCH ^oo;3015; 

ma> I MU- Q - T/C Press SAl 
MAX* THRUST 
coast 


pREL^'^’f^c^ ^0o;3O16; 

ma/T-u'^ q T/C Press SA2 

max* thrust 

Group Avg. % for: 

Prelaunch 
Max. Q 
Max. Thrust 
*Coast 


2A12C 

4 

/ 

• 0 

0 / 

.0 

0 / 

*c 

0 

/ 

.0 

A0120 

C 

/ 

.0 

0 r 

.0 

0 / 

.0 

0 

/ 

• ‘ f 

1P36C 

58 

/ 

.3 

•39 / 

.2 

31 / 

.? 

17 

/ 

• 1 

17 ^, 0 ? 

0 

/ 

*0 

0 / 

.0 

0 / 

.0 

C 

/ 

• c 

12060 

r \ 

/ 

• 0 

0 / 

.0 

0 / 

'•C 

7 

/ 

• 1 

200SC 

360 

/ 

1*8 • 

5 / 

.0 

0 / 

*0 

' 0 

/ 

♦ 0 

5680 

149 

/ 

1*5 

16 / 

.2 

7 / 

. 1 

2 

/ 

.0 

8800 

0 

/ 

• C 

0 / 

• 0 

0 / 

.0 

0 

/ 

.0 

12C60 

0 

/ 

*0 

0 / 

.0 

0 / 

*0 

4 

/ 

• 0 

E0060 . 

1615 

/ 

8*1 

546 -/ • 

2.7 

308 / 

1.5 

R2 

/ 

• 4 

9680 

366 

/ 

3*8 
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^ Table 7. Quiescent Signals only Active During Coast 
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Table 12. Percent of Transmitted Data 
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PERCENT SAMPLE 
NOT TRANSMITTED 

PROBABILITY OF 
STRING = n 

AVERAGE STRING 
OF ZEROS 

5 

100 

100.0 

5.00 

6 

97 

97.0 

5.97 

7 

88 

85.0 

6.85 

8 

80 

68.0 

7.63 

10 

73 

50.0 

8.82 

11 

68 

34.0 

9.56 

12 

59 

. 20.0 

10.05 

23 

57 

11.0 

11.47 

24 

97 

11.0 

12.85 

25 

88 

10.0 

14.06 

26 

80 

7.8 

15.00 

28 

73 

5.7 

15.74 

29 

68 

3.9 

16.25 

30 

59 

2.3 

16.57 

34 

57 

1.3 

16.80 

35 

99 

1.3 

17.03 

41 

99 

1.3 

17.34 

42 

97 

1.2 

17.64 

43 

88 

1.1 

17.92 

44 

80 

0.9 

18.15 

46 

73 

0.6 

18.32 

47 

68 

0.4 

18.44 

48 

59 

0.3 

18.53 

59 

57 

0.1 

18.57 

60 

97 

0.1 

18.61 

61 

88 

0.1 

18.65 

62 

80 

0.1 

18.69 

64 

73 

0.1 

18.74 

65 

68 

0.05 

18.77 

66 

59 

0.03 

18.78 

67 

57 

0.02 

18.79 

68 

98 

0.02 

18.80 

69 

97 

. 0.02 

18.81 

77 

96 

0.02 

18.82 

78 

97 

0.01 

18.83 

79 

88 

0. 01 

18.83 

80 

80‘ 

0.01 

18.84 

82 

73 

0. 01 

18.84 

83 

68 

0.005 

18.85 


Table 13. String Length Calculations 
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Figure 1. Signal Characteristics During Maximum Q 
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Figure 2. Average String During Maximum Q 
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APPENDIX B 
LEVEL I THEORY 


1.0 DESCRIPTION OF LEVEL I 
1*1 Level I Algorithm 

Level I encoding of binary digital data is applicable to amplitude-, frequency-, and phase - 
modulation systems as well as baseband systems. This appendix discusses the analysis 
of baseband systems employing Level I encoding* The methods utilized and the con- 
clusions drawn may be expected to be applicable with some modifications to modulation 
systems. 

In the baseband case, the information is encoded by causing transitions between two 
voltage levels at times determined by the encoding algorithm. To facilitate the discussion, 
consider the Level I encoder to put out an analog voltage v{t) that is either +V or -V. 

(See Figure 1. ) 

v(t) 



Figure 1, Level I Encoding 

The binary source produces a logic 1 or 0 every T seconds. (The information rate is 1/T 
bits per second. ) Now consider the time domain to be divided into cells (called bit cells) 
of length T seconds* The binary source associates a 1 or a 0 with each bit cell* The Level I 
encoder causes a transition between voltage levels ±V to occur at the beginning or at the 
middle of each bit cell, or no transition to occur, according to the following algorithm. 


I 

Bennett, W. R. & Davey, J.R, , Data Transmission, McGraw-Hill, N,Y. 1965, Chpt. 7-10 
and particularly Chpt. 11, 
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1. To encode the first bit cell (starting in either level with x representing either a 
0 or a 1)'J 

(i) If the first three source bits are llx or 000, encode a transition at the beginning 
of the bit cell. 

(ii) If the first three source bits are lOx, encode a transition at the middle of the 
bit cell. 

(iii) If the first three source bits are Olx or 001, encode no transition. 

2. To encode the bit cell (k s2): 

a. If the (k-1) bit cell had no transition, apply the rules l(i) through l(iii) to source bits 

k, k+1, k+2. 

st 

b. If the (k-1) bit cell had a transition at the middle of the cell, encode no transition. 

st 

c. If the (k-1) bit cell had a transition at the beginning of the cell, 

(i) If source bits k, k+1, k+2 are Ixx, encode no transition. 

(ii) If source bits k, k+1, k+2 are OOx, encode a transition at the middle of the 
bit cell. 


With one exception, these rules accomplish the encoding of transitions according to the 
information contained in pairs of source bits (with the appropriate constraints so that the 
information is encoded one-to-one). The exception is that three successive O's are made 
to cause a transition at the beginning of the bit cell and another transition at the middle of 
the next bit cell. This exceptional encoding is provided to assure transitions in the encoded 
bit stream when a long string of O^s is present. With this one exception, the rules cause a 
transition at the beginning of the bit cell if the source pair is 11, a transition at the middle 
of the bit cell is the source pair is 10, and no transition if the source pair is 01 or 00. 
Except for the three O's case, the constraint prevents transitions from occurring in 
successive bit cells. An example of a Level I encoded bit stream is given in Figure 2. 


SOURCE ^ 01 010000 1 110010 

data . 

LEVEL I BIT — —I— 1— 

STREAM I— 1— 1 . M I ,-J 

LOGICAL n- m- n+ m+ n- e- m+ n- e- n+ m+ n- n- m- 

CHARACTER^^^T -p 
IZATION ^ 
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Figure 2, Level I Encoding, Example 
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1«2 Level I Probabilities 

The Level I encoded bit stream will be analyzed as a stochastic process. The joint pro- 
babilities characterizing the process will be computed from the statistics of the binary 
source data and the Level I encoding algorithm. For this purpose, each encoded bit cell 
may be logically characterized as having an edge transition (e), a mid transition (m), or 
no transition (n) (as described in the algorithm), and, as having been in the upper level (+) 
or the lower level (-) at the beginning of the bit cell. (See Figure 2), (The analog signal v(t) 
making square transitions between voltage levels ±V, as in Figure 2, was introduced as a 
convenience in visualizing Level I and is not required to identify its statistical properties. ) 


The notation p(0j3, k) with a = e, m, or n, and = ±, will be used for the probability that the 
til 

k encoded bit cell is in the state For example, p(e+, 17) is the probability that the 
th 

17 bit cell has had encoded an edge transition beginning from the upper level. Similarily, 

the notation p( k; a k+1; ; k + s) is used for the joint probabilities, 

and k + n + l; ; , k + s | a jS*, k; k + 1 ; 

Q, (n) ^(n)^ k + n) fpr the conditional probabilities. The latter is the probability that the 


k + n + 1 


st 


th 

k + s^^ bit cells are in ths states a 


(h+l) ^(n+l) 


^(n + s)^(n+s) 


given that the k^^, , . . , k + n bit cells are known to be in the states oj , . . . , 
respectively. 


Throughout this analysis, it is assumed that the O’s and I's of the binary source occur with 

equal probability and that the information bits are statistically independent. In the case of 

statistically independent source data, the Level I algorithm makes it clear that the Level I 

th 

bit stream is a Markov process because the content of the k encoded bit cell depends only 

st 

on the content of the (k-1) . The conditional probability introduced above then satisfies the 
equation. 

|o^,k; a'jS’, k+1;. k + n) 
k+s k+ n) (1) 


p(Q^(^+l)^(*^‘'*‘l)^ k+ n+ 1;,..; k+s 

= k+n+l; 
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Finally, it is assumed that the starting time has receded to minus infinity so that the 
Level I bit stream may be analyzed as a stationary process* In that case p(a^) may be 
written for P(Q^ k), p (a/S, a ^ , . . . , jS for 

p,aAk, , k+ s), and p( I > 

for k + n + 1 ; k+s | k+ n). 

Stationary Markov processes are completely characterized by p(tt^ ) and the conditional 
probabiUty p{a’i3' |a/3). which is the probability that the bit cell in question is in the state 
a'j3' given that the previous bit cell was known to be in the state This fundamental 
conditional probabiUty \a6) may be thought of as the transition probabiUty from one 

bit cell to the next. To emphasize this interpretation, it is given the special notation 

p(a I ajS) = q(aj8 , O'’ ^ 

All the joint and conditional probabilities for the stationary Markov process may be written 
in terms of p(a /3 ) and q (o: ^3 , a’ j3 ' ). As an example, it may easily be shown that the 
order joint probability may be written 

( 3 ) 

The values of the transition probabiUty matrix elements may be determined by inspection 
using the Level I algorithm and the statistical independence and equi-probability of the 
source l"s and 0"s. They are 


e+ m+ n+ 

e+ pb 0 0 

m+ 0 0 0 

n+ 3/8 1/4 3/8 


e- 0 1/3 2/3 

m- 0 0 1 

n- 0 0 0 


* t 

e- m- n- *~0' p 
0 1/3 2/3"“ 

0 0 1 

0 0 0 


0 0 0 
0 0 0 
3/8 1/4 3/8 
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To find the values of p(a/S), it may be observed that p and q satisfy the equation 
L p{ajS) q (a d, a = p(a' 

2 

This equation may be obtained as a limiting form of the Smoluchowski equation. 

The values of the elements of p obtained on solving Eq. (5) are 

e+ 
m+ 
n+ 

p ( ) = e- 

m- 
n- 

1*3 The Probability of Sampling Next to a Transition 

In the system studied in Section 2 of this Appendix, the received Level I encoded wave train 
is first passed through a filter and then sampled in the middle of each half bit cell to de- 
termine which level the signal is in. The filter of Section 2 causes considerable rounding 
of the transitions (see Figure 2 where square transitions are shown) so. that the sample 
values (in the absence of noise) are ± 1/2 V rather than ±V if the sample is taken adjacent 
to a transition, it is, therefore, of interest to compute the probability that a given sample 
is taken next to a transition. 

Each half of the bit cell must be considered separately. For samples made in the middle of 
the first half of a bit cell, the information is contained in p(Oj(3), The sample is adjacent 


3/28 

3/28 

2/7 

3/28 

3/28 

2/7 


( 6 ) 


*M.C. Wang and G.E. Uhlenbeck, Rev. Mod. Phys. 1^, 323 (1945). Reprinted in 
Wax, "Selected Papers on Noise and Stochastic Processes," Dover, 1954. - Take the 
limit s “♦<* in Eq. (13) and use the discrete form of Eq. (4). (Note that Wang and 
Uhlenbeck use a sinistral convention in writing conditional probabilities, whereas I 
use a dextral convention except for q, ) 
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to a transition if there was an edge or a mid transition; otherwise it is not* Therefore 

Pr (adjacent) = p(e+) + p(m+) + p(e-) + p(m-) = 3/7 (7) 

l®^alf 

For samples made in the middle of the second half of a bit cell, p(Of^, Q^p ) is required. 
The sample is adjacent a transition if a mid transition occurred in the first bit cell 
(assuming the sample was taken in the first bit cell) or if an edge transition occurred in 
the second bit cell. The probability is, therefore 

Pr (adjacent) = p(m+, n-) + p(m-, n+) + p(n+, e+) + p(n-, e-) (8) 

2nd half 

r f 

Using Eqs. (3), (4), and (6), the values of p(a3* a /S ) are found to be 


e+ 

m+ 

n+ 

e- 

m- 

n- 

Thus 

Pr (adjacent) = 3/7 
2 nd half 


e+ 

m+ 

n+ 


m- 

n- 

0 

0 

0 1 
1 

0 

1/28 

1/14““ 

0 

0 

0 1 
1 

0 

0 

3/28 

3/28 

1/14 

3/2sl 

1 

0 

0 

0 

0 

1/28 

1 

i/i4 

0 

0 

0 

0 

0 

3/281 

0 

0 

0 

0 

0 

0 1 

3/28 

1/14 

3/28 


(9) 


( 10 ) 


Therefore, the result that will be used in Section 2 is that the sampler will measure 
±1/2 V (in the absence of noise) 3/7 of the time and ±V 4/7 of the time. 


2.0 ANALYSIS OF A LEVEL I SYSTEM 

2.1 Description of the Delco Laboratory System 

The system analyzed in this section corresponds closely to the system operating in the 
Delco laboratory. The system is shown in block diagram form in Figure 3. 
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V \ 



Figure 3. Delco Laboratory System — Idealized 

The main idealization in the block diagram is that the actual system contains additional 
elements for clock alignment and baseline stabilization between the receiving filter and 
the sampler. Their effect on performance is neglected in this analysis, and perfect clock 
alignment and baseline stabilization are assumed. 

The first three blocks of Figure 3 correspond to the two blocks of Figure 1. In Figure 3, 
the Level I encoder puts out a logical characterization, while the signal generator puts out 
the actual analog wave train such as is shown in Figure 2. The transmitting filter, line, 
and receiving filter cause a rounding of the square transitions shown in Figure 2, and it 
is the function of the sampler and threshold comparer to square the wave up so that it 
may be Level I decoded. The samples are taken in the center of each half-bit cell and 
compared against a zero threshold. 

When the line imposes a bandwidth restriction and it is desired to transmit information at 
the highest rate consistent with that restriction (which is the usual application of Level I), 
the information rate, iAt, is chosen to be the Nyquist rate; namely, twice the bandwidth 
of the line. The receiving filter is then chosen to have the same bandwidth (1/2T) so as to 
maximize the output signal-to-noise ratio. The transmission filter is usually omitted. 
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2,2 The System to be Analyzed — Notation 

To define the system to be analyzed, the signals generated together with the transfer func- 
tions of the various filters of Figure 3 are specified. The signals and transfer functions 
are not (necessarily) chosen to correspond to those actually used in the laboratory, but 
rather are chosen to optimize the output signal-to-noise ratio at fixed transmitted power 
and to minimize intersymbol interference. This approach allows more general conclusions 
to be drawn from the results than would otherwise be the case. The method followed is 
that of Sunde (which generalizes earlier work of Nyquist) as described in Bennett and Davey 
(Footnote 1, Chapter 7). 


For this analysis, the simplified block diagram of Figure 4 provides an adequate repre- 
sentation of the system. 



Figure 4. Simplified Block Diagram 


The block marked ’’Level I Source” combines the binary source and the Level I encoder of 
Figure 3. Its output is a sequence of logical symbols characterizing the Level I coded in- 
formation stream. For the present purposes, it is more convenient to use a different 
logical characterization than that shown in Figure 2 (that is, the characterization by 
a = e, m, or n; /3= + or -). Each half-bit cell (of length T/2)is, therefore, assigned a +1 
or a -1 depending on whether the voltage of Figure 2 is in the upper or the lower level 
during that half-bit cell. The +1 or -1 of the k half-bit cell is given the notation Sj^. 
Then the complete Level I encoded bit stream may be logically characterized by the 

sequence |s^l , k = -1,0, 1,2, The joint probabilities of a sequence |sj^| 

may be determined easily from the joint probabilities of Section 1. 

The block marked transmitter represents the signal generator, transmitting filter, and 
line of Figure 3. These functions may be lumped together without loss of generality. 

The transmitter is considered to put out a signal pulse or its negative in each half-bit cell 
depending on whether s^ is +1 or -1 in that half-bit cell. (The individual signal pulses 
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extend beyond the boundaries of the half-bit cells, of course, because the system is band- 
width limited). In setting up a notation for the signal pulses, it is convenient to relocate 
the time origin so that t = 0 occurs at the middle of the first half of a bit cell. (See Figure 
2 where t = 0 is at a bit cell boundary. ) Among other (notational) advantages, this reloca- 
tion allows the sample points to occur at t = k ^ , with k taking on integral values. The 
notation s(t) will be used for a (positive) pulse occurring in the zeroth half-bit cdll. The 
pulse occurring in the k^^ half-bit cell then becomes Sj^s(t ” k -^ )» and an expression for 
the whole transmitted wave, S(t), is 

<X> rj, 

S(t) = E s^s(t-k- ) (11) 

k=-oo 


To illustrate this notation, the data of Figure 2 are redrawn in Figure 5 for a typical s(t), 

T 

logically characterized by Sj^ and with s(t - k “^ ) shown in the zeroth and the eleventh 
bit cells. 



Figure 5. Typical Transmitted Wave (Figure 2 Data) 

Note that symmetrical pulses have been drawn. Following the treatment in Bennett and 
Davey (Reference 1), the analysis assumes symmetrical pulses s(t) = s(-t) and works with 
nonrealizable filters. This does not entail any loss of generality since the equivalent sys- 
tem function could be approached arbitrarily closely using realizable filters with 
sufficient delay 

The noise injected between the transmitter and the receiving filter is taken to be white 

with unilateral spectral density N . That is, if the noise voltage injected ds N(t), its mean 
^ o 

and autocorrelation functions are 
E [N(t) ] = 0 

N 

R (t-V) = E [N(t) N(t^)] = 6(t-t') (12) 
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The impulse response function of the receiving filter is called h(t), and the filter output 
for a given input is designated by affixing a subscript 1 to the input function* Thus, cor- 
responding to the single pulse input s(t), the wave train S(t) and the noise voltage N(t) are 

the outputs 00 

Si(t) = f dt' h(t - t’ ) S(t') 

-00 

00 

Sj^(t) =/dt' h(t - t') S(t’) (13) 

-GO 

00 

N^(t) =/dt’ h (t - t') N(t') 

-00 


It follows from Eq. (11) that 

00 ^ 
s. (t) = 2: (t - k — ) 

k=-oo 

and from Eq. (12) that 


(14) 


E [N^(t) ] = 0 

N 

R^(t - t') =E (N^(t) N^(t') ] = 


00 




dt' h(t - t") h (t^ - t") 


-00 


(15) 


To summarize all this notation, the block diagram of Figure 4 is redrawn and the functions 
defined listed on the appropriate part of Figure 6* 




Figure 6. Notation Summary 
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It is convenient in sonae stages of the analysis to work in the frequency domain. Let x(t) 
represent any of the functions of time defined above. The Fourier transform of x(t) is 
called X (to) and is defined by the relations 

x(co) = y dt e"^^ x<t) 

(16 

00 

= dco ^icot 

-J 2tT^ 

-00 


The Fourier transforms of Eq. (11) — (15) are 

coT 

_ 00 “ik 

S(co) = s (CO) S s e 
k=-(» K 

N 

Kto) = 

Sj^ (CO) = h(o) s (CO) 

S,(co) = h(c4 S (CO) 

,, coT 
00 -ik 

S,(a>) = i-tco) S s. e 
^ ^ k= -00 

Rl(co) = — I h(co) I 


( 11 ’) 

( 12 ’) 

(13’) 

(14’) 

(15’) 


R(C0) and R^(co) are, of course, the power spectral densities of N(t) and Nj^(t), respectively, 
h (co) is referred to as the receiving filter transfer function. 


2,3 Optimum Signal Pulses and Filter 

As mentioned before, the transmitted signal pulse s(t) and the receiving filter transfer 
function of the system to be analyzed are chosen to maximize the output signal-to-noise 
ratio at fixed transmitted power and to minimize intersymbol interference. Following 
Nyquist, the output single pulse spectrum s^(co) is chosen to be a raised cosine of band- 
width 1/2T. That is, 


- , , 1 2 coT 

Sj^(co) " 2 

0 




(17’) 
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or equivalently, 

Si(t) = i V 


sin 2TTt/T 
27Tt/T [1 -(2t/T) 


(17) 


This choice does not completely eliminate intersymbol interference with respect to samples 
taken in the middle of each half-bit cell (at t ~ k T/2); however, since 


s^(0) = I V ; (± I ) - I V ; 

Sj (k ^ = 0; k + 0, ± 1 


(18) 


it prevents intersymbol interference except between neighboring samples. For a Level I 
system such as the laboratory system in which (1) the received wave train is squared up 
by sampling in the middle of each half-bit cell and comparing this value against a threshold, 
and (2) information is sent along a bandlimited channel at the Nyquist rate, the intersymbol 
interference cannot be reduced further. 


Eq. (18) shows the meaning of the normalization constant V, If the nearest neighbors of 
a given half-bit cell both have pulses of the same sign, the nearest neighbor contributions 
add, and the sampled value is ±V* If the nearest neighbors of the given half-bit cell have 
pulses of the opposite sign, their contributions subtract leaving the bit cell's sampled value 
of ± 1/2 V. (See the comments at the beginning of Section 1,3.) This normalization is 
different from that su^ested by the labeling of Figure 2, which implies that ±V are the 
actual levels in the transmitted signal train. In the present convention, ± V are the levels 
(steady state after a long string of +1^ or -1®) after filtering by the receiving filter. This 
normalization is more convenient and, in any case. Section 2,5 shows how to e>q)ress V 
in terms of both the average and the peak transmitted power. 

With regard to the optimum receiving filter transfer function, Sunde has shown (Reference 1, 
Chapter 7) that, with respect to samples taken at the middle of the half-bit cells, the maxi- 
mum signal-to-noise ratio (with a white noise input) is obtained by choosing 

I h«o) I ~ I i (<o) I (19) 

The normalization constant can be chosen so that the filter has unity gain at zero frequency, 
and the phase characteristic can be chosen for convenience. Note that if the normalization 
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and phase were chosen so that h (a) = s (cu) the statement in the time domain would be 
h(t) = s(-t). Thus it is seen that Sunde’s result is just a generalization of the matched 
filter to the case where the signal pulses have non vanishing values outside their own half- 
bit cells. The optimum receiving filter (to within a nonessential normalization and phase 
characteristic) is, therefore, just a filter matched to the transmitted signal pulse function s(t)* 


From Eq# (13^) and (19) it may be inferred that 
s^ioo) - I h(co) f 

We therefore choose (see Eq. (17') ) 

«“ CoT I I 2 

h (CO) = cos — ^ ; 1 CtM < -ip 


Eq. (13'), (17'), and (19) then give 
s (CO) = |vT cos : I col 

0 ; I 6o| 


< 

> 


2jr 
T 
2 7T 
T 


or equivalently 



cos 2 77t/T 
[1 - (4t/T)^ ] 


( 20 ) 


( 21 ) 


( 22 ') 


( 22 ) 


This completes the specification of the system. To summarize, raised cosine spectra 
pulses were selected for the output pulses out of the receiving filter to minimize intersjrmbol 
interference, and the receiving filter was matched to the transmitted pulse shape to maxi- 
mize the sampled output signal-to-noise ratio. To list these results in a form for convenient 
reference. Figure 6 is redrawn and labeled accordingly in Figure 7, 


2.4 Sampling Theorems 


As indicated in Figure 7, the wave trains before and after filtering can be e3q)ressed in 
the forms 



(2S) 
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Figure 7. Section 2 'Level I System 


and 



(24) 


In this subsection, a number of sampling theorems are derived to enable these series to 
be summed, yielding simple analytic expressions for S(t) and S^^(t) whenever the series 
tOT 


k = -oo 


can be summed. The development of these sampling theorems proceeds as follows. 
First, the sequence represented by a discrete Fourier transform of minimum 


period T. gir/T 

s = / ^ 

J 277 


■hr ^ 

^ / 2 r 

e i (w) 


-277/T 


f - 2 S 8^ 

k= -00 


R75-34 


B-14 




\ 


DELCO ELECTRONICS DIVISION • SANTA BARBARA OPERATIONS • GENERAL MOTORS CORPORATION 

Then the function f(t) is defined according to 
217/T 

f(t) - / (26) 

- 2 ir/T 

T 

f(t) is simply a (convenient) analog function that goes through the points Sj^ at t = k - and 
contains no frequencies higher than I/'T, That is 

f ) = Sk (27) 


It may be constructed by associating the values s^ with the middle of each half-toit cell and 
then drawing the smoothest curve through the points* (See Figure 8* ) 



The usual sampling theorem may now be obtained by substituting Eq. (25) for f (to) into 
Eq. (26) and doing the co integration. The result is 


00 

f(t) = Z Sk 
k=-oo 


. 21T I T, 

sin — (t ^ k - ) 

f (t - k I ) 


(28) 


Eq* (2 8) is usually thou^t of as an analytical representation of a function f(t) that contains 
no frequency components higher than 1/T in terms of its sample values, where the samples 
are spaced T/2 seconds ^art. Here, it is regarded as the formula for an analog function 
containing no frequencies hi^er than l/T that has been constructed to have the property 


The next step is to express S(t) and Sj^(t) of Eq. (23) and (24) in terms of the function f(t) 
so constructed. To obtain a formula for S(t) in terms of f(t), note that (from Eq. (28)) 
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00 


4 ) = =^, 2 ^ \ 

^ Ir— -m ^ 


COS ^ <t - k “ ) 

T ^ 2 4 ' 


(29) 


After algebraic manipulation 


00 


f(t + 5 + t(t-T) = £ 


k=-oo 


®k 


2^^ - 1 , T , 

cos “Y ^ ^ 2 ^ 

TTp^‘] 


Comparing with Eq. (23) yields 


S(t) - |v (f(t + |)+ f (t -|) 1 


(30) 


(31) 


Similarly, to obtain a formula for S^(t) in terms of f(t), note that 

00 

2 


• 2ir .. , T , 

_ sin — (t -k- ) . 

f(t±?)=-E s. T 2 


k=-a, ±|) 

After somewhat more algebra, the result is 


\XJ 

|f(M-|) + f(t) + |f(t-|)= E 


Then comparing with Eq. (24) 


. 217 ,, T . 

sin Y” 2 ^ 


^(t-k|)[i-(|(t-k|)) j 




(32) 


(33) 


(34) 


The sampling theorems just derived, namely Eq. (28), (30), and (33), are not directly 
applicable to the analysis of Manchester biphase. To obtain sampling theorems that are, 
Eq, (25) is used to define a function T(t) according to 
417/T 

(35) 

zv ' ' 

-477/T 




which satisfies 


f’ (k 9 ) = s. 
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Thus f’(t) is an analytic function that goes through the points at t = k — , but it contains 
frequency components as hi^ as 2/T according to the prescription of Eq. (35 y used in 
conjunction with Eq. (25). The three sampling theorems that may be derived in terms of 
f'(t) are 


00 

f’(t)= L 3^ 
k= ~oo 


sin ^ (t - k ^ ) 
^ (t -kl ) 


(37) 


4^^ /f 1 T, 

rj, 00 . cos-;^ (t - k^ ) 

f(t^|)^f’(t-|) = 2: __Sk| 


8 


k= -oo 


■‘-(i.-'-if; 


477 X 

T rp 7 rp 00 sin Tj^t - k ^ ) 

if’ (t + i ) + f'(t) + i f'(t-i) = i 


k=-oo 4 




2. 5 Level I Wave Trains; Transmitted Power 


(38) 


(39) 


Valuable insights into the nature of the Level I link may be obtained by using the sampling 
theorems of Section 2.4 to get analytic expressions for S(t) and S^(t) in various special 
cases. The simplest case to analyze is that in which the original binary source produces 
all O's resulting in transitions in the Level I bit stream every 1-1/2 bits. For this case, 

3T 

f(t), S(t), ij, J dt S(t)^, and Sj^(t) 

0 

may then be computed. 


CASE 1 - Transitions every 1-1/2 bits. 



i 

► !• 
1 

1 

1 

1 

1 

1 

1 

• 1 

f 

1 

1 

1 

•1^ 

•1 

• 




1 

1 

!• 


1 

1 

I# 

1 

1 

1 

!• 



Figure 9, Case 1 Bit Stream 
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For this case. 


+ .... -l.O.l.l, 


(40) 


Therefore 


•1 ^ 

T ” - - 


f (CO) = 2 2 : s^e 

k=“CD 

,ujr 


(’ 1 ) 


k=-oo 


i3f^ i(3f+l)^ i(3f+2)^ 

e ^ +e ^ +e 


(41) 


T . ^2 icoxV y ^ f 

2 + e + e j L (- 1 ) e 


3coT 


The series S (-1) e may be summed according to the following method; 

f=-oo 


.. 3coT 

3T y ^ ft 2 
— H (- 1 ) e 

^ £=-<x) 


N-1 

lim, „ „ 2 

= N-tb ^ S e 

N odd ^ «=-(%^) 


U (-^ + ir ) 


lim. 
N-* 00 
N odd 


3T 

2 


.. .N-1. .3coT , 

-1 (-5~) ^ ) 


1 (n-1 )(2|^+ IT) 

1 + e +...+e 


lim. 

= N -+ 00 
N odd 


3T 




1-e _ 

i ^ aojl + 77 ) 

1-e 2 


(42) 


lim. 

= N -►00 
N odd 


2 


sin N ( + ir/2) 

sin ( + n/2) 


The diffraction grating amplitude function 3T/2 sin N (3co T/4 + tt/2) (x) [sin(3coT/4 + tt/ 2)] ^ 
has principal maxima at 3coT/4 + tt/ 2 = ntr, where n is an integer. In terms of (4 the 
maxima occur at Co= (2n-l) 2tr/3T.' In the limit as N goes to infinity, the function becomes 
a series of delta functions at these points. The normalization is determined by the normaliza- 
tion integral 
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(2n-l) ^ + 6 


lim. f 
N- 00 I 


2 


91T 

(2n-l) ^ 


sin + n/^ 

, /3C0 T 1T\ ~ 

^ 2 ) 


2Tf 


(43) 


Therefore, the result for the series is 


3T “ f 

“o A " (~f ) ® 

2 H= “OO' ' 


® r 2171 

27T L 6 CO- (2n-lh^ 
n= -00 -I 


(44) 


which, with Eq. (41), gives the result for the discrete Fourier transform of Sj^; namely 

. coT 

1 / 

f(co) = 


= |(l + e^ 2 +e‘“'^)2ir S dfw- (2 h-1) ||1 

\ / n= -CD ^ 


= 27T 


? i{ 

n = -00 \ 




)4' 


/o IV 2ir-i 
CO- (2n-l) ^ 


(45) 


Using the relation 


i(2n-l)tr/3 i(2n-l) 27r/3 

1 + e + e 


1 - e 


i(2n-l)7T 


1 - e 


i(2n-l)lr7^ ^ _ ^i(2n-l)7t/3 


(^6) 


Eq. (45) can be written in the form 

2/3 


00 

f(co) = 2iT L 


n= -00 1 - e 


1(2^1) ir/3 6 fo-(2n-l)f^] 


(47) 


The next step is the computation of f(t) using Eq. (26), The result is 

^ 277-1 . n\ 


f(t) = I 


sln4 


3T 

sin 7 t/6 


^ sin^^^+7r/2^ 


sin 7 t/2 


(48) 
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The factor multiplying the second term is 1/3 rather than 2/3 because the delta function 
has only been integrated halfway into. That this is the ri^t thing to do may be checked 
using summability procedures similar to those used in deriving Eq. (44). 


Finally, Eq. (31) and (33) may be used to determine S(t) and S^(t). The results are 
S(t) = I V cot 7T/6 sin^-^l^ ^ ^ 

Si(t) = I V cot 7T/6 cos 77/6 sin ^ \ (50) 


The average power in the wave train of Eq. (49) is 


3T 


3T 



dt S(t)^ 


2y2 

3 


(51) 


In cases 2-6 below, similar results are summarized for transitions every two bit cells 
through transitions every four bit cells (Figures 10 - 14). 


CASE 2 - Transitions every 2 bits. 


®k4 



i 

► !• 
t 

1 

t 

1 

1 

1 

1 


1 

1 

1 

1 

1 

k 




■ ■ ■— 

1 

1 

1 

1 



1 

1 

1 

1 





Figure 10. Case 2 Bit Stream 
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4T 



For Cases only the average transmitted power and Sj^(t) are listed, since S(t) and 
f(t) are rather simply related to S^(t) as in Cases 1 and 2* 

CASE 3 “ Transitions every 2-1/2 bits. 



Figure 11. Case 3 Bit Stream 



CASE 4 - Transitions every 3 bits. 



Figure 12. Case 4 Bit Stream 
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CASE 5 - Transition every 3--1/2 bits 



Figure 13* Case 5 Bit Stream 


7T 



CASE 6 - Transition every four bits 



Figure 14. Case 6 Bit Stream 
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On the basis of the above results for the transmitted power, it is clear that the average 

2 

transmitted power averaged over all possible Level 1 wave trains is between 2/3 V and 

7/8 V^- A straight average of these two values gives 0. 771 V^. The computer simu- 

2 

lation described in Section C computes the average power to be 0. 775 V . To eliminate 
lengthy decimal fractions from the results, the formula relating V to the average trans- 
mitted power is taken to be 

P = i (67) 

^AV 4 


This value differs from 0.775 by only 0.14 dB, which is negligible in the numerical 
calculations of Section D. 


The peak transmitted signal voltage required is that in Case 1 above. From Eq. (49), 
it is 2/3 V cot ir/l6 = V. The peak transmitted power required is therefore 


P = — 
^PK 3 ^ 



( 68 ) 


2.6 Level I Bit Error Rate 

(The characteristics of the system being considered are listed in Figure 7. ) Because of 

the structure of the raised cosine spectrum pulses, the sampler measures the sample 

values ± V when the sampled half-bit cell is not adjacent to a transition and ± 1/2 V 

when it is. As was shown in Section 1.3 , the samples are ±V 4/7 of the time and ±1/2 V 

3/7 of the time. The probability of making an error in comparing against a zero threshold 

to square up the filtered wave train is, therefore, 4/7 of the probability that the sampled 

noise voltage will exceed +V volts when the received signal voltage was +V, etc, , plus 

3/7 of the probability that the noise voltage will exceed 1/2 V, etc. The probability that 

2 

a Gaussian distributed noise voltage of zero mean and variance a will exceed a voltage 
level ±V is given by the standard formula 
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where 


erfc X = 


uu 

y-d^e- 


In the present case 


21T/T 


2 ^ r 6^ ^ o f <jtocos^-^= N /2T 

J 21T J 4 o 

-00 -21T/T 


The formula for the probability of making an error in determining the level of a given 


half-bit chosen at random is therefore 


Pr(c,4. 


■ 5 


Eq. (67) and (68) may be used to express this result in terms of either the average or 
the peak transmitted power. Thus 

Pr (€) =y • I erfc + f ‘ I 


4 1 




3 1 


p 

3 PK 


pr (e)= 3. iertc^f nW T M ^ 


The bit error rate (BER) is the probabilityof an informatioi bit being in error after Level I 
decoding. It may be obtained from Pr(€) by multiplying by 2 to convert from half bits to 
bits, and by multiplying by the number of decoded bits in error that result from one half- 
bit being in error (which will be called a). Thus 

BER = 2 O' Pr (e) (74) 

Arguments are presented below which suggest that 

a ^ 2 (75) 


(The actual value of oc is correctly taken into account in the simulation. ) 
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The complementary error functions in Eq. (72) or (73) are very rapidly decreasing func- 
tions of the arguments under the square roots. In all situations of interest, the first terms 
of Eq. (72) and (73) are completely negligible in comparison with the second terms. This 
implies that nearly all of the half-bit errors result from sampling the voltages ±l/2V; that 
is, in sampling half -bit cells that are adjacent to transitions. Thus, the half-bit errors 
almost always result in the shifting of a transition just one half-bit cell away from its 
correct position. Except in the cases where transitions are so close together that a shift 
in one results in a faulty interpretation of the other, a half— bit error that results in shifting 
a transition within the same bit cell (edge to mid, or mid to edge) results in one bit being 
in error in the decoded bit stream. A half -bit error that results in shifting the transition 
into a neighboring bit cell results in three bits in the decoded stream being in error. Thus 
on the average, neglecting the complicated cases, a single half-bit in error results in two 
bit errors in the Level I decoded bit stream. 

3.0 DIGITAL SIMULATION 

3. 1 Overall Simulation 

An overall block diagram of the digital simulation used to support this study is given in 
Figure 15. As illustrated in the figure, the simulation has been divided into three dis- 
tinct program steps. They are: 

X, Encoding Program 

2. Transmission, Receiving and Detection Program 

3. Decoding and Comparison Program. 

3.2 Encoding Program 

The purpose of the encoding program is to generate a magnetic tape containing the Level I 
coding of a random sequence of binary data. In addition, the output tape contains the 
original binary data. This data will be used to determine the number of errors made in 
the transmission, receiving, detection, and decoding process. 

A simplified flow diagram of the Encoding Program is given in Figure 16 A, and a FORTRAN 
listing of the program is given in Figure 16B, 
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Figure 15. Overall Simulation Block Diagram 
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Figure 16A. Program to Generate a Sequence of Random Binaiy Bits 
and Perform Level I Encoding Function 
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1 ; 

2'i 

"♦I 

51 

6 : 

7: 

8 l 

10 : 

Il- 

ls: 

13: 

U; 

15 : 

171 

18: 

19: 

c'O: 

21 : 

23: 

25: 
? 6 ; 
27: 
28 ; 
29; 
30: 
31': 
3 ^ : 
33 i 

36: 
’ 7 ‘l 
58; 
3g: 
40 : 
41 : 
42 « 
43: 
44 : 
45: 
46; 
47 1 
48: 
49: 
SC- 
SI: 
52: 
53; 
54: 
55: 
56; 
57- 
58! 
59: 
60; 
61; 
‘> 2 : 


REAL*4 mSq(15) 

DATA fiSG/' MID BIT TRANSITION N6 TRAkSTT^N '»• EDG5 TR 

lANSITION •/ 

INTEGER Data# BITCeUNt^ SHIFTCNT, DATATEMP, PITS/ DATABlTSiSEVEN 

INTEGER BIT0UT 

INTEGER 3IT0UTP 

DATA BITS/aZEOOOOOOO/ 

DATA SEVEN/7/ 

DIMENSI6N LEVELK64)/ LEVEL0UT(64) 

RECNT « 325 
IX* 84637 

CALL RANOeMI IX/ yPL ) 

IDATA « Ix 
BITC9UN”r • 1 
J ■ 1 

data » IDATA 

DATABITS • 1AND( ISA( IANO( DATA/ m TS )» -?^)/ SEVEN ) 
data - ISA( DATA# +1 ) i ShlFTCNT « 

2 IF{ DATABI TS*GE«6 .or* DATABITS*EQ- “• ) LFVtLI ( B I TCRUNT ) » +1 ; 

1 Gft TO 3 

IF{ DATABITS*GE*4 ) LEVELl (BITCGUNt ) s -1 ; r,6 T3 3 
LEVELl (BiTCeUNT) « 0 

3 DATABITS * IAND( ISA( IANCI DA^A/ B^TS U -?9)/ SEVEN ) 
data ■ TSAf DATA/ +1 ) t SHIFT^NT . SWIFTCf ^ + 1 

IF( 9HIFTcnT.GE« 29 ) GO TO llC 

4 BITCOUN'’' * BITCOUNT + 1 

IF( LEVELI (BITC5UNT-1) .E0«C ) Of TO 2 

IF( LEVELI (BITC6UNT-1) .n*-l ) I f VEL 1 ( B I TCOU^'T > = 0 ; SO TP 3 
6 IF( DATaBitS-GE» 4 ) LEVELl (BITCP'JNT) » 0 ; GO TP 3 
IF( 0ATA9 Its*LE. 1 ) LEVELl (-31 b : g*" TO 3 

WRI^Fi 6/200 ) ; ST 6 R 

POO FORMAT! TlO/ ’ AN UNACCOUNTFO FOR BIT pattfpN HAS OCCURFD') 
lln IF( ShIFTCNT .GT. 32 ) GO TB 121 
IF! SHIFTCNT .EQ* 32 ) GO TB 1 ?C 
IF! SHIFTCNT-GT.29 ) 60 TB 4 
100 CALI RANO0M{ IX/ YFL 1 
IDATaTMP a IX 

DATATE^^P - IAND( IDATATMP# BITS ) 

111 PATATEMP • IANO( ISA( DATATEMP/ -3 >, 8Z1C00O00G ) 
data « I0R( DATA/ OATATEMP J 
GB Tr 4 

120 BRI^E! 6^202 ) IDATA i WRITE! 7 ) IDATA 
IDAT^. « IDATATMP I DATA « IDATA 

G0 T- A 

121 wRlT^t 7 ) ( LEVELU I)/ 1*1/32) 

V WRI’‘E! 6/201 ) ( (MSG< (LEVELl ( I )-M )*5+K ) »Ksi/T)/ I»i,32) 

RECORD • RECSRD + 1 i SHIFTCNT « I ; BITCOU^:"^ » 1 
IF( record. EQ*RECNT ) END FILE 7 ; PPWIND 7 ; StOp 
CALL ANYPSONI 300S/ 1 ) 

301 CONTINUE 

CALL ANYPSONI 302S/3 ) 

IF! LEVEL1(32> .EQ. 0 ) 30 to ? 

IF! LEVELlOej.EQ.-l ) LEVELIU) » 0 | GO TO 3 
G6 T» 6 

300 OUTPUT! 102) REC0RD 
GB T* 301 

302 BUTPUT{102) RECARO 
ENDFILE 7 | REWIND 7 
5T0P 

?0l FBRMAT! T5/6C5A4)) 

20P FBRMAT! TlO/ 'RANDBM DATA INPUT STRE AM » / 5X. Z«/ ) 

E ND 

Figure 16B, FORTRAN Listing of the Encoding Program 
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3, 3 Transmission, Receiving and Detection Program 

This program represents the main thrust of the simulation effort. As indicated by the 
title, this program simulates the transmission process, receiving process, and detection 
procedure. The resultant output is a magnetic tape containing the detected Level I coded 
pulse train and the original binary data which it represents. 

While the transmission filter is explicitly defined in terms of its impulse response, the 
receiving filter is implicitly defined by the impulse of the "overall'^ system (transmission 
filter - receiving filter). 

The first step in simulating the various system outputs is to convert the input transition 
sequence to a discrete pulse train of four pulses per bit cell, i.e. , leading edge, one- 
quarter, half, and three-quarter points of the bit cell. This discrete pulse train is then 
convolved with the impulse response representing the transmission filter. This output 
is used to define the average transmitted power. 

The output of the overall system (receiving filter output) is defined by convolving the dis- 
crete pulse train with the impulse response of the overall system. As part of the simula- 
tion of the overall system output, it is necessary to account for the effect of noise. This 
is accomplished by generating band-limited white noise with zero mean and variance 
(as prescribed in Sections 2.2 and 2.6) and adding it to the already generated system 
output. This forms the overall system output to be used by the various detection schemes. 

The output representing the one-quarter and three-quarter points of each bit cell are now 
used as input to the detection procedures. Two procedures are suggested in Sections 2.1 
and 4. 3 and a simplified description follows; 

1. Single Point. This procedure issued a ±1 for each point depending on 
whether the amplitude is greater than or less than zero. The polarity 
of these two samples are now compared with the last previously gener- 
ated data and a decision made as to whether an edge, mid bit, or no 
transition has occurred. 
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2. Four Point. The second procedure forms two weighted values utilizing 
the present samples and the three immediately previous samples. The 
absolute value of these weighted values are then compared to a threshold. 

For those values which exceed the threshold, the local maximum is de- 
termined, defining the point at which a transition has occurred. 

An output tape containing both the decoded transition sequence and the original binary 
data is now generated. This output tape will be processed by another program to determine 
how many errors were made. A simplified block diagram of this program is given in Figure 
17A and its corresporfding FORTRAN listing is shown in Figure 17B. 

3.4 Decoding and Comparison Program 

This program takes the transition sequence generated by the transmission, receiving, 
and detection program and regenerates a binary data stream. This resulting data stream 
is then compared with the original binary data in order to determine how many errors 
were made as a result of the simulated system. 

Figure 18A is a simplified block diagram of this program, and Figure 18B contains the 
FORTRAN listing. 
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Figure 17 A. 
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PROGRAM 

INITIALIZATION 
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GENERATE THE TIME IMPULSE RESPONSE 
OF BOTH THE TRANSMISSION FILTER 
AND THE OVERALL SYSTEM 


READ IN ORIGINAL BINARY DATA AND 
THE CORRESPONDING TRANSITION SEQUENCE 

i 

TRANSFORM THE TRANSITION SEQUENCE 
INTO DISCRETE PULSE TRAIN WITH 4 POINTS 
PER BINARY BIT 

"X" 

COMPUTE OUTPUT OF TRANSMISSION FILTER 
TO CALCULATE AVERAGE TRANSMITTED 
POWER 

1 

COMPUTE THE OUTPUT OF THE 
OVERALL SYSTEM 

1 , 

GENERATE BAND LIMITED WHITE NOISE 
AND ADD TO SYSTEM OUTPUT 


SAMPLE THE MIDDLE OF EACH HALF BIT 
AND RECONSTRUCT THE LEVEL 1 
DATA STREAM 

1 

GENERATE AN OUTPUT TAPE CONTAINING BOTH 
THE ORIGINAL BINARY DATA AND 
LEVEL 1 CODED DATA 


Simulation of the Transmission, Receivir^ and Detection Process 
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u 

2 - 

3 ! 

4: 

5: 

6r 

7 : 

8 : 

iOj 

111 

12: 

13 ; 

14 ; 

15 ; 

16 : 

17 : 

18 : 

19 * 

20 : 

21; 

22l 

23 : 

24 : 

25 ] 

? 6 : 

27 : 

28 : 

29 : 

30 : 

31 ! 

33 : 

34 ! 

35 * 

36 : 

37; 

38! 

391 

40 ; 

41 - 

42 ; 

43 ; 

44 ! 

4S5 

‘•41 

47 : 

48'i 

49; 

50 : 

51 : 

S2: 

53 ; 

54 : 

55 - 

56 ; 

57* 

58; 

59: 


integer DATA/ BITSNE/ DATABITS# R 

DI^5EMSI9N X{41)> 0ATABUT(12)/ H(41), P(2>/ DA(6)i SIG{6)/ 

1 DATAeUTl ( 12)/ SCRIPTSI41), CaPS(41)/ Q(2) 

’ / L'^VELlN{32)/ 0K(2) 

dimension LEVELK128)/ LEVELPUPS)/ lEVEL3fl?8)# lEVEL4(64)^ 
♦ LEV‘^L5(64)/ LEVEL6(64) 

DATA BIT5NE/8Z80000000/ 
namelist SIGMA/ THRESHEj thresmdf 
IN6ISE ■ 49323 

IX ■ 84637 ; R » 9 ; BITCNT • 0 
DIFf^P • 0.0 / uEVEL3P *1*0 
LEVC-L2P * 1«0 i DIFFPl *0.0 
THRFSH2 * 2*5 i THRESHDF • 1-25 
IStGN'2 ■ 1 
D8 60 I - 1/41 

60 SCRTPTS( I ) • 0*0 
SCRIRTS(21) * 4/3*1415927 
SCRIPTS (20) * SCRIPTS (22) < 1.0 
DB 61 I » 2/20*2 

SCRIPTS( 2l- I ) f 4*(-*l )■*#(!/?)/( 1-I**2}/3«1^15927 

61 SC»IPTS(2l+I ) « SCRlPTS(2l*I ) 

DC 62 I ■ 1/41 

62 H( I ) » 0*0 

H(21 > • 1.0 

H(?C) s H(22) . 5/(3.1415927*3) 
hi(3.9) • H(a3) p 0*5 
DO 63 I * 3/19/2 

H(?:i-I) * (-l)**(I/2* - 1/2. ) *8/(4-T**2)/3. 1415927/3 

63 H(21+I ) ■ H(21-l ) 

I Input 

eUTRl.'T SiGMA/ THRESH2/ THPES^DF 

« 0.0 

CALL STRIP'JLSE(20) 

CALL STPlPSPEED( 10) 

D8 30 I • 1/6 

30 DA( n * 1.0 

CALL STRlPL5T(DA) 

DO 31 I » 1/6 

31 DA(I) ■ -1*0 
CALL STPIPLOT(DA) 

D6 32 I • 1/6 

3? DA( n » 0.0 

CALL STRlPULSEdBO) 

CALL STRlPL8T(OA) 

CALL STRIPULSE(IO) 

PLAST • 1.0 i SIGNLAST = 1*0 
D0 10 V ' 1/325 

READf 7, END-lOO ) IDATA ; RFAD( 7) (LEVELIN( J) * J»l/32> 

D6 9 J » 1/32 

IF( LEVELIM( J) .LT.C ) p(l)» PLAST i P(2) *-p(D 
IF( LEVELIN(U) *EQ.O ) R(1>:;* P(2) * PLAST 
IF( LEVELIN( J)*GT.G) pm » P(2) • -PUAST 
PLAST ■ P(2) 

DAU) » ,LEVELIN(J)/2.0 

1 DU 2 M • -20/20/1 

2 X(M+?1) • X(M4-21) + P( 1 )tH(M^21 ) 

DATAfiUT(R) • XU) j DAT^fltUT(R4l) * X(2) 

D© 4P MM * -20/20/1 

42 CAPSIMM+ai) * CAPS(MM+21) * P( 1 )*SCRIPT$ ) 


Figure 17B (Sheet 1 of 4). FORTRAN Listing of the Transmission, Receiving 
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60 ! 
61 ! 
62: 
63; 
64: 
65; 

66j 

67: 

68 : 

69; 

70; 

71; 

72; 

73: 

74: 

75* 

771 

78* 

79: 

80 : 

81; 

82; 

83; 

84: 
85i 
861 
87: 
88 : 
89* 
90 ; 
91: 
92; 
93; 
94: 
95: 
96* 
97: 
9a; 


4,3 


P(2)#H(M+21) 
i DATAOUTCP+S) 


X(2) 


44 


45 


P8i^T- - CapS( 1)**2 + CAPS(2)##2 ♦ P9WER 
Nip B IMP + 2 ' 

DO 3 K. » -20>18#1 
X(K + 21) *■ X(K+23) 

■ X(41) # 0*0 
D9 43 KK a -20M8#1 
CAPSfKK + 21) • CAPS(Ki<-»-231 
CAP5(40) « CAPS(41) s 0.0 
D6 4 M » -20.20/1 • 

X(P+21> « X(M + 2D + 

DATA6UT(P+2) ■ X(l> 

DB 44 MM * -20/20/1 

CAPS(MM + ?i) k CAPS(MM-»-21) + P(2)#SCRIPTS(MM+?1 ) 

PBPEP " Cai^S(1>**2 + CAPS(2)**2 + P9WER 

NP » NP 2 

DB 5 K ■ -20/ 1 o/ 1 

X(K>?1) ■ y(K423) 

DB 45 Kv* > -20/18/1 
CAPSfKK+2l) = CAPS(KK+23) 

CAPSt4C) « CAPS(41) a C.O 
X(4C) » X(41) a 0.0 
DA(2) • P{l)/2*0 
DAO) • OataCjUT(R)/2.0 

CALL GAUSseM( IN0I6E. SIGMA/ XMEAN, SNPtSE 1 
DATAeUT(R) + SN0ISE 


OATARUTl(R) = 
call GA!jSS0M< 
DATA'>uT 1 (P+1 ) 
CALL GAUSSOMt 
DATAPUTl (R+2) 
CALL GAUSSBM( 
DATA^Ut’i (R*3) 


XMEANi 


SNBISE 
+ SNdlSE 
XMEAN/ 9N9ISE 
$N6ISE 

XMEAN, SN9ISE 
+ SN0I5E 


> 


) 


) 


I NOISE/ sigma, 
a DATA9UT(P+1> 

IN8ISE/ SIGMA/ 

« DATAPUT(P+2) 

IN0ISE/ sigma, 

« 0ATAPUT{RO) 

DA(A) > DATAeUTl(R)/2.0 
IF( BITCNT.LT. 5 > OB T9 51 

IF(ABS(DATABUT1(R) ).GT.0.UOK(naSIGN(l.O/DATAdjTl(R) ) 

IF( aBS(CataBUT1(R+2) ) .GT.0.1 )GK(2).SIGN(1.0/0ATABUT1(R^2) ) 
LEVELBUT • 0 

IF($tgnlAST.EQ.OK( 1 ) .AN0.SIGNLAST.E3.0X(?) ) 
IF(SIGNLAST.NE.QK(1) .AND.SIGNLAST.NF,0K{2) ) 
IF{STG''LAST.EQ.QX( 1 ) *ANCi,$IGNLAST.NE*flK(2) ) 


LEVEleUT 
LEVEL0JT 
LEV flout 


0 

1 

-1 


99;' 

SIGNLAST - 0K(2) 


100; 

LEVEL4(3lTCNT-4) - LEVElOUT 

101: 

0A(5) « LEVELBUT/2.0 


102: 

IFt abs(SIGO)-OA(5) ) .GT 

. 0*1 ) DA(6) • 0.5 1 ERRBRCNT 

103: 

* ERRORCNT -k 1 


104: 

IF(APS(0ATABUT1 (R) ) .GT.O. 

1 )LEVEL1 ( {BITCNT-5)*2+l ) = 

105* 

* SIGN!a.O,DATA-iUTl (R) ) 


106; 

IF( ABS( DAT ABUT! (R+2) ) .QT* 

0.1) LEVELl ( (BITCNT-5)*2+2) * 

107; 

# SIGN(1.0,DATAeUTl(R+2) ) 


108: 

DIFFI a 0.0 ; DIFF2 a o.r 


109': 

Dirri » AbS( OATABUTUP) 

- DATA0UT1 :R-2) ) 

110; 

DIFF? a ABS( DATAeUTl(R+2 

) - DATAOUTl(R) ) 

111: 

if: DIFFl *LE* THRt.St^DF ) 

DIFFl a 0*0 

112; 

IF( 0IFF2 .LE. THRESHDF ) 

0IFF2 a 0.0 

113: 

IF( DIFFPl .QT. THRESHDF 

) G6 TP 120 

114* 

IF( DIFFl .GT. THRESHDF ) 

50 T*’ l2l 

iiSi 

122 LEVEL2M8ITCNT-5)*2+1 ) = 

LEVEL2P 

116: 

IF( DIFF2 .GT. THRESHDF ) 

GP T0 124 

117; 

123 LLVEL2( t?^ITCNT*5)#2 + 2) a 

LFVEL2P 

118: 

DIFFPl ■ 0.0 : 


119‘: 

G8 T0 126 



A 

Figure 17B (Sheet 2 of 4). FORTRAN Listing of the Transmission, Receiving 
and Detection Process 
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1?0! 

i?l: 

122: 

123] 

124 * 

125: 

126 ; 

127: 

12S': 

129^ 

130; 

131; 

13?; 

133; 

134; 

135; 

136; 

137; 

138; 

139; 

140J, 

141: 

142; 

143; 

144; 

145; 

146; 

147; 

148; 

149; 

150: 

151: 

152: 

153 ; 

154; 

155; 

156; 

157; 

158; 

159: 

i<^0; 

161; 

162: 

163] 

164; 

165: 

166 ; 

167; 

168: 

169] 

170: 

I7ii 

172] 

174; 

175: 

176] 

177: 

178 ] 

179: 


12*^ IF( OlPTPl .GT* Dim ) r,e Tfi 125 
LEVEL2( ('^ITCNT-5)*2) » LEVEL2P 

121 IF< Dim tGT. ABS( DIPPPI) .AND* DIFFPl .LT ,0* 0 ) GO T0 129 
IF( Oim *LT* DirF2 ) GO T6 122 
LEVEt2nBlTCNT«5)*2 + l) - -LEVEIPP 
LEVEL2P “ -LEVEL2P 
LEvEL2( (BITCNT-5>*?+2) » LEVELPP 
OIFFPI » -DIFF2 
GO TB 126 

129 LEvEU 2( fBlTCNT-5)*2-l) • -LEVEL? ( ( B I TCNT-5 ) *2- U 
GO TO l?2 

125 IF( DIEP2 ,0T. DIFFt #AND* 0TFF2 .GT. THRESHOF ) G6 T8 12g 
LEVEL2 ( (BitCNT-5)*2) --LEVEL2P 

LEVEL2P " -LEVEL2P 

LEVEL2( (BitCNT-5)*?+1 ) ® LEVEL2P 

G6 TO 1?3 

128 LEVEL2( (8ITCNT-5)*2) » LEVEL2P 

LEVEL2( (BitCNT-5)*?+1 ) - -LEVEL2P 

LEVEL2r- » -LEVEL2P 

LEVEL2( (5ITCNT-5>#2+2) « LEVELPP 

DlProi . 0.0 
GB TO 126 
124 DlFPPl « DIFF? 

126 CONTINUE 

iFf LEVFL2( {BITCNT-5 )*c + 1) .NE* ISIGN2) LEVE1.2BT « 1 ; Gt' TQ j ?7 
IF( lEVEL2( <BITCNT-5)*2>1) .EQ*LEVEL2( (HTTCnT-5)*2+?) )LFVFL2?t.o ; 

* GO l27 

LEVE120T « -1 

127 IStO‘i2 » LEVEL2( {BITCNT-5)*2+2l 

IF< ABS{ SIG(2) • LEVEL28T/2.0) .GT.0.1 ) ERpCNTI • EPRCf^TI 1 
0IFF3 e 0IFF4 s 0*C 

DIFF9 p AbS( DATA3UTURJ + 0 .5*0 ATA^UT 1 ( R-? > - 0*5*DATA9I1T1 ( P-A ) 

* - DATA8UTKP-6) ) 

r'IFF4 » ABS( DaTA 0UT1(P^2) + 0.5*DATA9UT1IP) - D.5*DATA»"T1 fP-2) 

* - DATA9UT(R«4) ) 

IF{ DIFFP *GT. THRESM2 ) G0 TB 110 
IF{ DIFF3 *GT. THRf.SH? ) GO TO 111 
11? LEvEL3( (B1TCNT-5)*2*1) . LEVEL3P 
TF{ DIFF4 .GT. THRCSH2 ) GO TO U4 

113 LEVEl3( (BITCNT-5)*2+2) • LEVEL3P 
OIFFP * 0.0 

Gfi 116 

no IF( DIFFP .GT. DIFF3 ) GO TO 115 
LEVEL3( (BITCNT-5)*2> » LEVtL3P 
111 TF( DIFF3 .LT. DIFF4 ) GO TO 112 
LEVEU3( (BITCNT-5)*2+D » -LEVEL3P 
LEVEL3P • -LEVEL3P 
GO TO 113 

115 LEvEL3( (BITCNT-5)*2) * -LEVEL3P 

LEVEL3P * -LEVEL3P 
LEVEL3MBITCNT-5)*2 + 1 ) -LEVEL3P 
GO TO Il3 ^ ■ 

114 DIFFD • 0IFF4 

116 CONTINUE 

IF( pITCNT,lT.6 > 05 TO 51 

IF{ LEVFL3(tBITCNT. 5)*2).NE. LEVELS! (8ITCNT,5>*2-1 ULEVELGOT • 1 t 
1 GO T8 117 

IF( LEVEL3( (8ITCNT-$)#2).EQ. LEVELS! (PITCNT*5>*2+1) )LEV£L30T • 0 ; 
1 GO TO 117 
LEVEL30T • -1 
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180: 117 IF( ABS( SIGU) • LEVEL3RT/2*)fGT*0»l) EPRCK'T? • ERRCNT? f 1 

181] LEVEC6t6lTCNT-5) • LEVEL36T 

182J 51 SIGH ) « SIG<2) 

1835 SIQf?) B S1G(3) 

1845 SIG(3) ■ SIG(4) 

185] S1G(4) • SIG(5> 

186; SIG(5) « 9IGC6) 

187; SIG( 6 ) » 0A(1) 

188; y call STdTploT(DA) 

189| DAO) « DatA 6UT<R-»-1 )/2i0 

190; DA(4) * nATAeUTl{R4l)/2,0 

19.1J X CALL STRlPLRKOA) 

192'; DAO) . P(2)/2.0 

193; DAO) * DATAGUT(R+2)/2tO 

194; DA(4) • DATAeUTl(R+2)/2*0 

195; X CALL STPlPLeT(DA) 

196; DAO) B DATAeuT(R+3)/2*0 

197; DA(4) B nATA8UTl(R+3)/2*0 

198; X CALL StpIpU9T(DA) 

199; BITC^'T * BITCNT + I 

200; oe 11 I » lf 8 

SOI; DATAOUT(I) • OATAdUT(W^) 

202; 11 DATARUTKI) • DATAeUTl ( 144 ) 

203; DA( 6 ) ■ 0.0 • 

204; 9 CeNTiNUF 

205; if; ^ -LT. 2 ) GB to 131 

206; lr’F<rTE( 8 ) ( LEVEL4(L)^L« 1>32 ) 

207; WRITE(9) { LEVEL 6 (L)/L* 102 ) 

208; D 6 130 KK . 1|32 

209; LEvE|^4(KK) «. LEVEL4(KK*32) 

210; LEVEL4(KK432) 0*0 

211; LEVEL 6 (Ki^) • LEVEL6( KK+32 ) 

212; 130 LEVEl6(KK+32) • 0-0 

213; oe 132 KK « 1,64 

214; LEvELSfl^^) • LEVELS! KK+64 ) 

215; LEVELS * 0*0 

216; LEVELSfXK) • LEVEL3 t KK+64 ) 

217; 13? levels ( xX4b4) • 0*0 

218; BITCMT B BITCNT - 32 

219; WBRDCNT * WBRDCNT 4 1 

220; I3l WRITE!?) IDATA 
221; WRlT|(q) idATA 

222; 10 ceNTINUE 

223; 100 08 3? I " 1^6 

224; 33 DA! n » 0,0 

225; X CALL STPlPULSE!200) 

226; X CALL ST«= IPL0T! DA ) 

227 ; PBWERAVO • P 8 WER/NP 

228; BITCN'T « BITCNT ♦ W0RDCNT*3c 

229; BUTPgT pBWERAVG, ERR 6 RCNT# ERRCNTl, BITCNT 

230; 8 UTPUT PRRCMT2 

231; REWIND 7 

332; WRITE! r) I LEVEL^!L)# L ■ 1<32 ) 

333; WR1TE!0) ( LEVEL 6 (L>/ L ■ 1*32 ) 

234: ENDEILE 8 J REMIND 6 

235; ENDFILF ? J REWIND 9 

236; STBP 

237; END 
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1: BitCM 

2: DIMENSTf'N tEVELlIN<64)# LEv'EL9UT^64) 

PEaD( 7/END«100 ) lOATAl i RFAD( 7 ) ( LEVELl IK( J)# 3^) 

REaD( 7>EnD*100 ) IDATAS t REAOt 7 ) ( LEVf L 1 1 Ni ( J ) j J«33# A4 ) 

5: Nf « 1 

6: IDATACL'T « 0 

7: 3'* Db 1] * N/3? 

8: IF ( k-.r ".I ) Gb T9 32 

9: 33 IF( LEVrLlIM(K-l >.GT*-1 ) GO T9 30 

10: IF( LEvFLIIM(K) .GT-0 ) WPITFI A, 300 ) 

11; LEVEUe'^'T (K) ■ 0 i G0 TP 40 

12: 300 F6FMtT{ Tlo, t A MIO eiT IRaNSITIR^ MpT F'^LLOi/fO BY A ?ve ^RANFITjR^f) 

13: 3r> 1F( I^EYFLI IN’(K-1 ) •E3«0 > of' T9 31 

1^: IF( UEv'f^Ll TN'(k) ,FO.O ) LFVtL'^O'T (K ) « 1 ; T6 40 

15: IF( LEvFLI INtK) .EQ*-1 ) LEVFL9'JT(K) « 0 ; 30 T9 

16: WRITE! A# 301 ) i Lt.VELbUT(<) * 0 ; 39 Tf» 4r 

17; 301 FfcR^ATr Tjq, tywO CS'^iSECUT I vF EDGE TRANSITIONS OfTfCTEH') 

18; 31 IF( IFvrt- 1 If- (K-1 ) ♦N'E'O ) WRITE! <-t302 > ; IFVEL'^UTfO « 0 j 

is: * Of T9 40 

20: 3? IF( LEvrL.l ri(K) •E0«0 ) LE L 9UT ( k ) » 0 ; G‘ T0 40 

21: 1F( l^E'vTLlINtK) .EQ.-l ) LEvELBUKK) » 1 ; GO TO 4c 

22: IF( UEvrLl IN(K+1 ) *EG«C ) LEVEL91 =T(k1 • 1 ; G9 T5 4C 

23: IF ( LEvTLI INIK+l) .EQ»-1 ) LEVELPUT(k) « 0 ; GR t0 i^q 

24; WRITC( a» 301 ) j LEV^LouT(^) » 0 ; SP T9 4r 

25_; 302 FOPMa"^ ( TlO» ' AN ILLEGAL BIT tPA,n'SIT!‘’K C90E gftEcTEC) 

26: 4r IF( IFLAC.fq.I ) GO TO 52 

27: OATATFMf « ISC( LEVELPUT(K),32-k ) 

28j ll IDaTACUT = ISRI IDATAOLT# DATATr*^P ) 

29: DATATEmo . IE9R( IDATAl^ I0aTA9UT ) 

30; IF( 0ATA’T«P-E3*C ) GO T9 5C 

31; 08 51 M 8 1/3? 

32; ERRQRCWT a IAN0( DATATFMP/1 ) + FRR^RCNT 

33; 51 OATATE''!^ - ISC( OATATENP/-1 ) 

34; 50 EITCNT » SITCK’T + 3? 

35; K » 33 ; Iflag« 1 

36; G6 TR 33 

37: 5? IFLAG > 0 ; N ■ 2 

3S: 0ATATE^^c , isc( LEVEL9UT ( 33 ) i 3t ) 

39; IDaTAOUT * 0 

40; IDATaP-JT = I0R< 1DATA8UT/ DATATEMP ) 

41 J Oe 53 J = 1/32 

4?: 53 LEVELlIf-U) » LEVELIIN! j+32) 

43; IDATAl 8 IDATA2 

44: REAO( 7*END8101 ) IDATA2 ; READ ( 7/ ENO ■ 101 ) ( LEVEL 1 IN ( J ) /j«33/ 64 ) 

45; G6 T9 34 

46; 100 WRITE! 6^103 ) i ST9P 

47; 103 F8RMAT! 10/ 'END 9F FILE ENG9UNTERED AT BEGI^'MINg gp TAPE ' ) 

48: 101 WRITp! 6/104 ) BITCNT / ERR9RCNT ; ST9P 

49: 104 F8Rf"AT{ Tio/ls/ i BITS WPPE PROCESSED '/ 1 5/ ' ERRORS ^^ERE DPTECTED'l 

50: END 


Figure 18B. FORTRAN Listing of the Decoding and Comparison Program 




R75-34 


B-37 



OELCO ELECTRONICS DIVISION • SANTA BARBARA OPERATIONS • GENERAL MOTORS CORPORATION 


4.0 COMPARISONS AND DISCUSSION 

4.1 Results for Manchester Bi-Phase 

A Manchester bi-phase system may be defined and then optimized using the Nyquist-Sunde 
approach just as was done for Level I. In Section 4.1 the results (without any derivations) 
are given for a Manchester system with error detection capability. Such a system sacri- 
fices some margin over noise in order to detect (without correcting) the majority of errors 
that occur. For this purpose, the signal train is sampled in the middle of each half -bit 
cell and compared against a zero threshold to determine the signal level in each half-bit 
independently. The intent is to have single half-bit errors show up as nonvalid Manchester 
pulses. 

The optimum Nyquist-Sunde system is shown in Figure 19 . The notation used is virtually 
identical with that of Figure 7. The Manchester encoding algorithm is simply that s^= +1, 

= -1 (k even) if the corresponding binary source bit was a logic 1, andSj^ = -1, 

= +1 if the binary source bit was a logic 0, 



Figure 19. Manchester Bi-Phase System 
The series for S(t) and Sj^(t) were summed in the following two special cases, 
CASE 1 - All Manchester 0- Phases 
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Figure 20, Case 1 Bit Stream 
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~ (-1) ; k — “1,0, 1,2,. •• 

S(t) = /^ V cos ^ 

c. /4-V T7 2^ 

S^(t) = V cos 


(76) 

(77) 

(78) 


CASE 2 - Alternate Manchester 0- and 77“ Phases 
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Figure 21. Case 2 Bit Stream 


^ 211—1 ^ 2 £ ^ •••» “ 1 , 0 , 1 , 2 , ... 


S(t) = /2 V 


S^(t) = V 


cos 77/8 cos 


77/^ + cos “l^cos 

1^^^) 


The average transmitted power in both Case 1 and Case 2 is 

P = 

^AV ^ 


(79) 

(80) 
( 81 ) 


(82) 


and the simulation confirms that this is the correct identification. The peak power re- 
quired is that which occurs in Case 1 and is 


P 


PK 


2v2 = 2 


(83) 


Finally, the bit error rate is 
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4,2 Manchester Without Error Detection 

If the error detection capability is relinquished, a more than 3 dB improvement in margin 
over noise may be effected for Manchester bi -phase. The simplest method of obtaining 
this improvement is to compare the difference between the first half-bit sample and the 
second half-bit sample against a zero threshold to determine whether the full bit was a 
zero phase or a pi-phase The effective signal power is thereby quadrupled while the 
noise power in the difference is only doubled. The effective signal-to-noise ratio is thus 
increased by a factor of two. The probability of error does not have to be doubled to 
convert from half bits to bits, and so the bit error rate is 



Eq. (85) may be recognized as the bit error rate of an optimum performance communica- 
tions system operating at bit rate l/T in the presence of white Gaussian noise of spectral 
density The Manchester system of Section 4,2 is, therefore, an optiinum system. 


4,3 



The improvement obtained on using the difference of Manchester half-bit samples as a 

test statistic suggests that Level I systems employing sample differences might obtain 

an improvement over the system of Section 2. The simulation has been used to carry out 

a (nonexhaustive) study of such differencing schemes. The most successful of the schemes 

studied employed a linear combination of four neighboring half-bit samples to determine 

when the transitions between the levels take place. If the voltage output of the receiving 

T 

filter, S, (t) + N. (t), is called V. (t), then the sampled values are V^(k - ). (See Figure 7.) 

XX X • X ^ 

The statistic computed to determine whether a transition occurred between the k-1 and 
the k^^ half bit cell is 


Gk - I Vi ((k^l) I) ^ I (k I ) - I ((k-1) I) - ((k-2) I ^ I 


As the Level I wave train undergoes a transition, successively takes on the values 0, 
i V, 7/4 V, 5/2V, 7/4 v, V, 0 (in the absence of noise). The 5/2 value is identified as the 
transition point. Identifying the transition point in the presence of noise is treated as a 
combined detection and estimation problem, Gj^ is compared against a threshold of 3/2 V 
in order to detect that a transition has taken place. To estimate the location of the 
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transition points, the values of Gj^ above the threshold are examined to identify the local 
maxima. The results of the simulation are discussed below in Section 4*4. 


Another Level I system that was studied for the sake of comparison was the system of 

Figure 19, with the Manchester source replaced by a Level I source. Since in this case, 

2 

all the sample values are ± V (and ~ V ), the bit error rate is (compare Eq. (72), 


(74), and (84) ). 
BER = 20f 



(871 


This system will be referred to as the "no intersymbol interference" (No I-I) Level I 
system in the discussion below. 


Alpha is again the number of bits in error in the decoded Level I bit stream per half-bit 
error developed in the threshold comparison. The computation of a is now, however, 
more complex than in the arguments following Eq. (75). When a half-bit error occurs 
adjacent a transition, we may again ejq)ect approximately two bit errors to develop for 
each half-bit in error on the average. Now, however, since all the sample values are 
±V, a significant number of half-bit errors also occur not adjacent to a transition. These 
errors generally result in nonvalid Level I wave trains and must be eliminated before (or 
during) Level I decoding. The Level I decoding algorithm may therefore be considered to 
correct many of the half-bit errors that occur not adjacent to a transition (the number 
corrected will be a function of the algorithm), and a will be significantly less than 2. 
Some crude estimates of O' were made using the joint probabilities of Section 1 and a 
simple assumption about which errors are corrected. The results suggest that 

I £ a < 1 (88) 

The simulation results for No I-I Level I presented below include a simple algorithm 
for obtaining a valid Level I wave train before decoding, and therefore include effects such 
as possible dependence of a on signal- to-noise ratio. 


Still another Level I system was identified as worthy of study, but its analysis could not 
be completed within the scope of the current contract. It is shown in Figure 22. 
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Figure 22 • Level I , Transition Detection 


The output pulse function 

Sj^(t) =^os 2irt/T |i -^t/':i^'j 
has the property that 

(±T/4) = 1, ^ (k T/2 - T/4) = 0 for k ^ 0, 1. 

Since the samples here are taken at the edges rather than at the centers of the half-bit 
cells, this system has the property that the sample values (in the absence of noise) 



th til 

are ±V if there is no transition between the k-1 and the k half-bit cell and 0 if there 
Is a transition* One way of detecting transitions, then, would be to compute the test 
statistic 

where V^(t) is again S^{t) + Nj^(t) and compare against a threshold of 1/2 V. 
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Finally, it may be noted that block and sequential decoding schemes have not been studied 
even though they might be expected to provide significant improvement in Level I per- 
formance because of the correlation present in the Level I bit stream. Such studies would, 
of course, be well outside the scope of the current contract. If or when they could be 
carried out, they would probably identify systems with better error performance at the 
cost of increased hardware complexity. Such systems would vitiate some of Level I's 
principal advantages; namely, simplicity and low cost. 

4.4 Comparisons 

The bit error rate of the Level I system studied in Section 2 is plotted in Figure 23, to- 
gether with the BER's of the Section 4 systems (Manchester and Level I), The range lO"*^ 

-5 

to 10 bit errors per bit is shown. The argument is taken to be 10 log Pat/ZN /T, which 

is the ratio in dB of the average signal power to the noise power in a bandwidth equal to 

the information rate. The curves based on the analytic expressions derived are drawn 

solid while the curves drawn through points computed with the simulation are dashed. The 

simulation results were computed using, at most, 10,368 information bits, so that the 

lower points on the simulation curves have a greater variance. ^ 

The curve labeled ''Manchester Without Error Detection" is a plot of Eq, (85), As men- 
tioned in Section 4.2, it is an optimum performance curve. No system can obtain better 
performance (in white noise) without going to block or convolutional coding techniques. 

The solid curve labeled "Manchester Bi -phase" is a plot of Eq. (84), The dashed curve 
was obtained with a Delco simulation. Although the development of this Manchester simu- 
lation was not part of the contract, the techniques employed were very similar to those 
of the Level I simulations. The curve is .therefore shown to demonstrate the close agree- 
ment (better than 0.1 dB) between the simulation and the theory. 

The solid curve labeled "Level I (Section 2)" is a plot of Eq. (74) using Eq. (72) with 
a = 2. Recall that a is the average number of bit errors that result in the Level I decoded 
bit stream due to a half-bit error. (See the arguments following Eq. (75) in the text. ) The 
simulation curve could be fit pretty closely by Eq. (74) with a= 1.4. In fact, an exact 
fit could be obtained by an a that varied continuously from 1.2 at the upper end of the 
curve to 1. 5 at the lower. 
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In comparing the curves, it is seen that a Level I, Section 2 system has between 1 and 
1-1/2 dB less margin over noise that the Manchester bi-phase system of Section 4*1 * 
To comment in detail and summarize the specific effects that lead to this result, refer 
to Eq, (84) for Manchester, namely 


®E^Man = 2- 


and to Eq* (74) with the first term of Eq, (72) neglected. 



(84) 


I 

The first terra in Eq, (72) has been neglected because the ±v samples don’t contribute 
significantly to the error rate over the range considered. The 3/7 factor is present be- 
cause the ±V samples only occur 3/7 of the time. The factor gl is discussed above and In 
the text following Eq, (75), Referring to the arguments of the square roots, there are 
three effects which result in Manchester having a 3/2 (or a 10 log 3/2 = 1. 8 dB) advantage 
over Level I in signal-to-noise ratio. 



• The Level I filter has half the Manchester filter bandwidth, so the 
Level I detection process is only contaminated by half the noise. 

• The Level I samples are ±1/2 V compared with ±V for Manchester; 
therefore, on a relative basis. Level I has 1/4 the signal power. 

• Signal-to-noise ratio comparisons are made on the basis of average 

transmitted power. For Level I, ~ 3/4V^ while for Manchester 
2 A V 

= V . Therefore, in converting to an average power basis for comparing 
signal-to-noise ratio, Manchester loses a factor 3/4 with respect to Level I 


Putting these three effects together, Manchester is seen to have a 1/2 x 4 X 3/4 = 3/2 
advantage over Level I in signal-to-noise ratio. 

The points on the curve labeled ’’Level I (four-point difference)" were computed with 
the simulation of the system described in Section 4.3. The four-point difference scheme is 
seen to have equal or poorer performance than the single-point Level I system of Section 2. 
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This result can be understood qualitatively as follows. The statistic of Eq. (86) has three 
times the noise variance of the single sample. In estimating the effective signal power, 
the detection and estimation part of the problem must be treated separately. In detecting 
a 5/2 V value by comparing against a threshold of 3/2 V, the effective signal power is , 
an improvement over single-point Level I by a factor of 4. In estimating the location of 
the maxima, however, we are normally comparing 5/2 V with 7/4V. One way to do this is 
to compute the difference and compare against a zero threshold. In computing the difference 
the noise variance is again doubled, while the effective signal power is now 

(5/2V - 7/4V)^ = 9/16 V^. 

A crude estimate of the effective signal- to-noise ratio compared to single -point Level I 
is therefore 

1/3 X 4 X 1/2 X 9/16 = 3/8, or 10 log 3/8 - 4. 3 dB. 

The arguments are admittedly crude and, in fact, give a considerably poorer estimate of 
four-point differencing than is actually achieved on the simulation. On the basis of the simu- 
lation results, and to a lesser extent on these arguments, it is concluded that simple 
differencing schemes do not allow a significant improvement in Level I error performance. 

The curve labeled "Level I (No I-I)" is drawn through the points generated by the simu- 
ulation for the "no intersymbol interference" Level I system discussed in Section 4. 3 . 

Eq. (87) gives an anal3d:ical result of the bit error rate for this system. This equation 
was not plotted in Figure 23 becaise for a = 1, it is identical with the Manchester bi-phase 
curve. Eq. (88), and the arguments preceding it suggest that 1/2 s a ^ 1, and the curves for 
other values of a can easily be inferred from the curve for a = 1. Figure 23 shows two 
of the simulation points lying near the Manchester curve and the others lying above it. This 
suggests a rather strong dependence of Ofon signal- to-noise ratio. (Although the behavior 
of a, and therefore the Level I error correcting capability, depends on the particular de- 
coding algorithm chosen, a study of decoding algorithms and attempts to optimize a were 
considered beyond the scope of the contract.) In any case, the Level I (No I-I) per- 
formance was, at most, 0. 8 dB poorer than Manchester Bi-phase. It is felt that with proper 
algorithm design, this type of Level I system could be made to perform at least as well as 
Manchester over a specified range of BER. This is an improtant conclusion. Although 
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the Level I system of Section 2 is optimum for its usual application in which it is desired 
to send information at the highest rate consistent with the bandwidth restriction, the No 
I-I Level 1 system of Section 4, 3 obtains a 1 to 1-1/2 dB Improvement In margin over 
noise by standing the bandwidth by a factor of two. In applications in which this band- 
width expansion is allowed, the following statement may be made. The bit error rate 
performance of Level I is comparable to Manchester Bi-phase (with error detection) and 
only 3 dB poorer than Manchester without error detection. 
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APPENDIX C 

EVALUATION OF DELCO ELECTRONICS 
3-PHASE DPM AND CONVENTIONAL 4-PHASE DPM 

1.0 INTRODUCTION 

This appendix develops an expression for the phase error probability as a function of 
signal- to-noise ratio for L-phase differential phase modulation and presents an evalua- 
tion of Delco Electronics 3-phase DPM (differential phase modvilation) technique versus 
conventional 4-phase DPM. 

The results of the phase error probability analysis are shown in Figure C-1. As can be 
seen, at a si gnal-to- noise ratio of 10 dB, the probability of a phase error in 3-phase 
modulation is down by a factor of 10 from that of 4-phase. This factor increases with 
increasing SNR, and the advantages of 3-phase become more prominent. Translating 
phase error rate into bit error rate shows; (1) for conventional 4-phase, a single phase 
error produces, on the average, one bit error, while (2) for Delco's 3-phase, a phase 
error produces, on the average, approximately one and one-half bit errors. Comparing 
the two schemes at 10 dB, then, the 3-phase bit error rate is down by a factor of about 
seven from that of 4-phase. 

2.0 DESCRIPTION OF THE MODULATION SCHEMES 

A binary stream is used to modulate the phase of a carrier in a differential manner. 
Figure C-2 shows conventional 4-phase DPM, where 00 causes a phase increment of 45^, 
01 a phase increment of 135^, etc. Figure C-3 shows Delco’s 3-phase DPM, where two 
consecutive I's cause a phase increment of 180®, two consecutive O's, a phase increment 
of 300®, and 10 a phase Increment of 60®. The 01 combination is encoded by holding the 
previous phase for one bit time, thus coding the initial 0 bit, with the 1 being coded with 
the next bit in the data stream. An example of coding with conventional 4-phase and with 
Delco’s 3-phase is shown in Figure C-4. 



PROBABILITY OF PHASE DIFFERENCE ERROR 
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Figure C-2. Conventional 4-Phase DPM Figure 3. Delco's 3- Phase DPM 


225° 45ojl35°jl35^45°j3lS°jl35«^15«j315<^5°|45°|225°|3150^^^ 

1 1000101 00100110 1 0 1 00011 10 1 BINARY STREAM 

|l80°j 300° |60°| 60° j 60° 180° 1 60° 1 60° |300°|l80° 60°| * -^ DE LCD'S 3-PHASE 

Figure C-4, An Example of Encoding in Conventional 4-Phase and 
in Deleons 3- Phase 
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3.0 PHASE ERROR PROBABILITY 

This evaluation is made by first presenting a derivation of the probability of a phase dif- 
ference error when the phase change between two signaling intervals is preserved as a 
reference (differential phase modulation). The traditional assumptions are made; namely, 
of an ideal digital phase detector observing a sine wave with one of L-phases in additive 
gaussian noise. 

The probability density of the error in phase, 0, is given by; 

2 

p(e)= 2 ^ e"^ [^i+yrp cose 6^“°® ® (i+erf(/p cose))] -IT (i) 

where p is the signal-to-noise ratio and el‘f denotes the error function. A derivation of 
Equation 1 is given in Section 5 of this appendix for completeness and because, even 
though it is reportedly a well known result, it is incorrect in the reference. 


Errors are committed when the detector measures a phase outside of the region from 
0 - i^L to 0 + 7t/L, where 0 is the received phase in the absence of noise. 


The single error probability, P(e), in DPM is thus given by the probability that the 
absolute value of a phase difference exceeds 7T/L| that is, 

P(e) = Prob (r/L < | - Sq I ) 

In order to derive an expression for P(e), it is assumed that the probability of $ lying 
between and 6^^ + dQ^^, and of 0^ lying outside the interval from B^-tt/Ij to 
0- + 7 t/L, is the product; 


p(e^) <3% 


,+ff/L 


P(e„) de„ 


0j^-ir/L 


J 


Letting 6 vary from -it to +JT and recognizing the symmetry of p(6) about 0=0, yields 

0 + ir/L T 


ir 

P(e) = 2 fp(6j) 
o 


eJ-Tt/h 

O 


P<9o> 


d0j^ 
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Since phase differential is used to compute the transmitted phase in DPM, successive 
errors are not independent, that is 


^(^2 I / 0 


for finite S/N, From the above, it follows that the joint probability of successive errors 
is given by 


7T 




r 01 +ir/L 

r 

1-/ P(V^o 


l-y P(92)<»2 



. 0j-n/L 




or in more compact notation 


TT 


P(e^,e2)= 2 Jp(0) 


0 +n/l. 

1 - I p(a) da 
0 - 1 T/h 


-/ 


d 0 


( 2 ) 


Recognizing that the conditional probability P(C2 j e^) is given by the ratio of the joint to 
the single error probability, we have 

P<"2lV= p|e,) 

and, thus, expressions for all pertinent probabilities. 


We now show that the conditional error probability of Equation 2 reduces to the results 

2 

of Salz and Saltzberg for L = 2, namely, 

n 


P(e^,62) =1 J 1 - erf (v/p cos 0) I p<0) d0 
o 




(3) 


Letting 0 + ^/2 

q(f^) = J p(a) da 
0 -- 13 / 2 

and differentiating with respect to 0 
dq _ 


^=p(0 + n/2) - p{0 - ir/2) 


R75-34 


C-5 



DELCO ELECTRONICS DIVISION • SANTA BARBARA OPERATIONS • GENERAL MOTORS CORPORATION 

Since 

cos (0± tt/ 2 ) = ^ sin 
we can write 



By making a change of variable 
X = ^ cos 0 

and recognizing that q(7r/2) == 1/2, we have 


cos 0 


q(^^) = 


J 


e ^ dx 


-00 


Since the integrand is an even function and 
00 


Tij 


e ^ dx = 1 


-oo 


Therefore, 

1 - q(0) = 


00 




dx = ^ [1 - erf i/p cos 0 ) ] 


/jo cos 0 

which proves the equality of Equation 2 and Equation 3 for L = 2. 


Figure C-1 is a plot of the single error probability for L - 2, 3, and 4. The evaluation 
of the Integral 

01 + tt/l 

f P(6^) 

- 7T/L 

requires numerical integration for the signal-to-noise ratios of Interest. The double 
integration was carried out using Simpson's rule and an error function approximation 
accurate to nine digits. 
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4.0 BIT ;errqr probability 

To evaluate the effectiveness of Delco Electronics 3-phase DPM versus conventional 
4-phase DPM, the process by which phase errors are converted to bit errors must be 
examined. In conventional 4-phase (Figure C-2), considering nonadjacent phase 
vector errors as highly improbable, phase errors translate to bit errors one-to-one. 
Three-phase modulation as shown in Figure C-3 results in a higher bit-to-phase error 
\ rate. Considering first the no-hold combinations, a "1 1" could be interpreted as either 
a’ ’TO” or ”00” with equal probability, which results in an average of one and one-half 
bit errors per phase error. The same is true for a 00 Interpreted as a 11 or 10. A 10 
Interpreted as a 11 or 00 results in one bit error per phase error. When a 01 follows a 
11 or 1 0, the phase angle is held for one additional bit time. To investigate the bit- 
to— phase error rate for this case, the six combinations of five binary digits beginning 
1101, 0001, or 1001 must be analyzed. There are 18 ways that these holding operations 
can result in bit errors. Considering these, it can be determined that the 18 possible 
phase errors produce 33 bit errors. Holding operations for random data occur with a 
relative frequency of one-fourth. The bit-to-phase error ratio for Delco’ s 3-phase is 
thus 



Thus the probability of bit error for Delco 's 3-phase is a factor of about one and one- 
half higher than the probability of phase error shown in Figure C-1. For biphase and 
4-phase, the probabilities of phase and bit error are equal. Comparing 3-phase and 4- 
phase at 10 dB then, the 3-phase bit error rate is down by a factor of about seven from 
that of 4-phase. 


5,0 DERIVATION OF PROBABILITY DENSITY QF PHASE IN PHASE MODULATION 

The received wave is assumed to have the form 
V(t) = (A + x) cos (ojt + 0 ) + y sin (cot + 0 ) 

where 0 is the phase that would be measured by the ideal phase detector at the end of 
a t-second signalling interval in the absence of noise, and x = x(t), and y = y(t) are the 
inphase and quadrature Gaussian noise components, respectively, each with average 
noise power of CT^. The phase of the received wave Is then 

Phase of V(t) = 0 + tan“^ 
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The next step is to find the probability density function (pdf) on the error in phase, p(0 ), 


where 


e = tan’^ JL- 

^ A+x 


Letting 

u = A + X 

V = y 

the pdf on u and on v are given by 

-(u-A)^/2o-^ 

f(u) = e 


and 


JW a 

- 2 2 
1 -V /2a 
g(v) = — e 

/W a 


Recognizing that 8 = tan"^ v/u, one can transform to polar coordinates, (u - r cos 6 
and V = r sin 6) to find the joint density function in r and 8, q(r,8 ), with p(8) then being 


given by the Integral 
00 

P(0) = J^ q(r,6) dr 


Following this procedure, we obtain an expression for q(r, 0 ) 


-A^/2a r -'(r'^“2rA cos 9 )/2a^ 


q(r,8) = 


2 ir a 


re 


] 


Letting the signal-to-noise ratio be 


p - A^/2a^ 


and 


2 1 ,_.2 ^ .2 2 


(r - 2rA cos 8 + A cos 0 ) 


2a 


the pdf on phase error is given by 


00 


pO) = 


1 g -P (l-cos 8) 


/2ira 


f -s2 

/ ( /%a s + A cos 8) e ds 

-y/p cos 8 
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which becomes an integration 

p(6) = |l + /Tp cos ® (1 + erf ( /pcos B ) )J 

where erf denotes the error hinction 

a 

erf (a ) = — 

/¥ 
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r -X dx. 

J ® 
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APPENDIX D 

CONVOLUTIONAL ENCODING AND VITERBI DECODING 

1.0 INTRODUCTION 

This appendix presents; (1) the essential features of convolutional encoding and Viterbi 
decoding, (2) a listing of a digital computer program used to evaluate the compatibility 
of Level I encoding and decoding with the above mentioned techniques, and (3) the results 
of simulations with random data sequences. 

2.0 CONVOLUTIONAL ENCODER 

A general binary convolutional encoder consists of a bK stage binary shift register and 
V mod-2 adders. Each of the mod-2 adders is connected to certain of the shift register 
stages. The pattern of connections specifies the code. Information bits are shifted 
into the encoder shift register b bits at a time. After each b bit shift, the outputs of 
the mod-2 adders are sampled sequentially yielding the code symbols. 

An example of a simple convolutional encoder with K = 3, b = 1, and v = 2 is shown 
below. The binary input data bits 0, 1, 1, 0, 1, 0... generate the code outputs 00, 11, 
01 , 00 , 10 , ... 


010001 



011010 .... 

DATA SEQUENCE 


4 

V 


R75-34 


D-1 






DELCO electronics division • SANTA BARBARA OPERATIONS • GENERAL MOTORS CORPORATION 

The state of the convolutional encoder is the contents of the first b(k-l) shift register 
stages. The encoder state together with the next b input bits uniquely specify the v 
output symbols* 


In the above example, the encoder state is specified by the contents of the first two shift 
register stages. Denoting the four possible states as a = 00, b = 01, c = 10, and d = 11, 
the input bits 0, 1, 1, 0, 1, 0... then correspond to the successive encoder states a, b, 
d, c, b, c * . * 

The allowable transitions between states at k - 1 and k are shown in the diagram below 
together with the transition code sequence. In this example, two paths lead to each of 
the four states. In general, there are 2*^ paths entering each of the 2°^^ states. 


a= 


c= I 10 I 



10 


Assuming a binary symmetric channel (i* e. , symbol errors are independent and occur 
with probability (p) and that all input data sequences are equally likely, a maximum 
likelLhood decoder is one which examines the error-corrupted received sequence and 
chooses the data sequence corresponding to the transmitted code sequence that is closest 
to the received sequence in the sense of Hamming distance; that is, the transmitted sequence 
which differs from the received sequence in the minimum number of symbols. 

A brute-force maximum likelihood decoder would calculate the Hamming distance on all 
paths and select the minimum. The information bits corresponding to that path would 
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form the decoder output. For an L bit information sequence there are 2^ paths in the above 
example. This method obviously quickly becomes impractical as L increases, 

3.0 VITERBI DECODING 

The Viterbi decoding algorithm greatly reduces the effort required for maximum likeli- 
hood decoding by taking advantage of the fact that the minimum distance path to, for 
example, state a at time k in the above encoder example can be only one of two candidates: 
(1) the minimum distance path to state a at time k-1, and (2) the minimum distance path 
to state c at time k-1. A comparison is performed by adding the new distance accumulated 
in the k^^ transition by each of these paths to this minimum distances at time k-1. 

In general, a Viterbi decoder calculates the likelihood of each of the 2^ paths entering 
a given state and eliminates from further consideration all but the most likely path that 
leads to that state. This is done for each of the states; after each decoding opera- 

tion, only one path remains leading to each state. In cases where distances are identical, 
an arbitrary selection is made. It is important to recognize that in eliminating the less 
likely paths entering each state, the Viterbi decoder wUl not reject any path that would 
have been selected by the brute-force maximum likelihood decoder. 

The great advantage of the Viterbi decoder is that the number of decoder operations per- 
formed in decoding L bits in only which is linear in L. The Viterbi decoding 

technique is limited to relatively short constraint length codes (K^ 6) due to the exponen- 
tial dependence of decoder operations per bit decoded on K, 

It has been shown that, with hi^ probability, the decoder- selected paths have a 

common stem, which eventually branches off to the various states. This further suggests 
that, if the decoder stores enou^ of the past information bit history of each of the 
2 b(K-l) oldest bits on all paths will be identical. It has been demonstrated 

theoretically and through simulation that a value of the length of the information bit path 
history per state of four or five times the code constraint length, K, is sufficient for 
negligible degradation from optimum decoder performance. 
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4.0 SIMULATION PROGRAM 

A digital simulation program was developed incorporating rate 1/2 convolutional and 
Level 1 encoding of random source data, error injection capability. Level 1 and Viterbl 
decoding, A listing of this program follows at the end of this appendix. 

5.0 CQMPAT ABILITY OF CONVOLUTIONAL AND LEVEL I ENCODING 

A number of computer simulations have been made in an attempt to determine the com- 
patibility of convolutional and Level I encoding. Two programs were employed. The first 
permits rate one-half convolutional encoding of random source data, error Injection, and 
Viterbi decoding. The second, discussed above, incorporates rate one-half convolutional 
encoding of random source data followed by Lievel I encoding, error injection, Level I 
decoding, and Viterbi decoding. The Level I decode uses the essential features of the 
Viterbi algorithm. Bit memory paths of length 16 and 32 bits were used in the Level I and 
Viterbi decoding, respectively. 

Each run was made with a source data sequence of 4,500 random bits. Error rates were 
established which were high enough to result in errors in the decoded sequence. Based on 
the models of commiinlcation links discussed in Section 4 of Appendix B, it was concluded 
that the probability of a channel error in the case of rate onehalf convolutional transmission 
is given by 

Pr (e) = ^ erfc V^N^ 

whereas in the case of rate one-half convolutional plus Level I transmission, the prob- 
ability is determined by 

Pr(e) = I 

From these equations, a signal-to-noise ratio which results in 700 errors in 4,500 bits 
of source data for rate one-half convolutional transmission while results in 1,420 errors 
in rate one-half convolutional plus Level I transmission. 

The results of simulations to test the compatibility of the two encoding schemes is shown 
in Table 1. For the lower error rates used, it can be concluded that the convolutional en- 
coding and Viterbi decoding is more effective in error removal. At the higher error rates, 
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the results tend to indicate that the inclusion of Level I encode-decode does not affect 
the error correcting capability of Viterbi decoding. 


RATE 1/2 CONVOLUTIONAL ENCODER 
VITERBI DECODER 

RATE 1/2 CONVOLUTIONAL AND 
LEVEL I ENCODER 
LEVEL I AND VITERBI DECODER 

— ■ ■ 

CHANNEL 

ERRORS 

CHANNEL 

ERRORS 

ERRORS 

AFTER DECODE 

ERRORS 

AFTER DECODE 

700 

73 

1,420 

308 

800 

129 

1,540 

400 

900 

194 

1,640 

439 

1,000 

204 

1,760 

459 

1,100 

382 

1,840 

512 

1,200 

421 

1,940 

532 

1,300 

516 

2,040 

564 

1,400 

585 

2,120 

560 


Table 1. Results of Simulations for Compatibility Test 
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SIMULATION PROGRAM 


l: 

21 

3 : 

4 ; 

5 : 

6 : 

7 : 

' 

lo; 

U; 

1 ?: 

13 ; 

14 ; 

1 * 5 : 

16 : 

17 ; 

lA* 

19 : 

20 : 

2 xi 

23 ; 

23 ; 

24 : 

25 : 

26 ; 

27 ; 

28 ; 

29 ; 

30 ; 

3l; 

32 '; 

33 * 

34 : 

35 ; 

36 : 

37 : 

38 ; 

39 ; 


c THTS PR0GRAM \S F0R The: C0NV9LUTI 8NAL ENCeOlNG OF 

C binary digits LISTCD in the matrix R> where K«3/ 

C B*l> Vp2 are the CSNVBLUTiaMAL PARAMETERS 

C 

IMPLICIT InTEGER(A-Z) 

C0MM9N Al (640)# A2(640)#PS(640)iDR(640)#Re(640) 
dimension raw(400.) 
dimension R(160) 

integer R 

REaD(5#100) (R{ I)j I«1/160) 

100 FeRMAT(? 0 li) 

call CENC9[)(R) 

c 

c Pack ai#ap inT6 raw 

D5 5 I«1/3P0#2 


RAWn )» Al( 1/2 + 1) 

5 RAw(I+l)» A2(I/2+l) 
call encode (320#RAW) 

C EPR0R iNTR^DuCTieN 

C ERR0R INTR9DUCTIBN 

integer ERR0R 


IX«67 

WRiTE(6/172) 

172’ ESRMATdHl) 

D0 1002 I-l#n 


call RAND9m< ix#yfl ) 

ERR9R »m9D( IABS(IX)#250) 
WRITE(6j 1003) ERReR#lX 
1003 F8RMAT(* Al A2 ERR9R AT» j2l4) 

IF(M0D( IXj 2) *£0.0) Al (ERR0R 

IF(MP0( lX,2)tEQ.l) A2 ( error 

1002 continue 


CAllt DECQDE (320#RAW) 

G SPREAD RAW BACK TO A1,A2 
DO 6 dl#320i2 

AKI/2+1 ) » RO (I) 

6 A2(I/2+l ) ■ RO (I + D 

call VDECOD(R) 

END 


) »N0D( Al (ERROR) +1/2) 
)«M0 D( a2(ERR0R) +1/2) 


0 / 

0 ^ 
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i: 

?; 

3: 

m : 

• 

6 ; 

7: 

8; 

9 : 

10 : 

li: 

1 ?: 

13; 

U: 

155 

16 ': 

17; 

18: 

19; 

20 ; 

21 ; 

29 : 

23 ; 

24 * 

25: 

2 $; 


C 

c 


200 

^00 

300 


subroutine CENC0D(R) 

IMPLICIT INTEGERIA^Z) 

CQMMBN A1 (640)/ A2(640)/PS(640)/DR{640)iRe(640) 
dimension R(160) 
integer Al 
integer A2 

THE SEQUENCE 81(1)82(1) REPRESENTS THE CODE FOR R(n 

81(1) » R(l) 

A2(l) • R(l) 

81(2) « R(l) + R(2) 

82(2) f R(2) 

IF(81(2) »EQ*2) A1(2) * 0 
D0 200 I’liil58 

AKU2) • R(I) + R(I + 1) + R(I+2) 

A2( T+2) “ R( I ) + R( 1+2) 

IF (81 ( I) *£0*2) 81 ( I) « 0 

IF( 81 ( I ) *EQ#3) Aid) * 1 

IF(A2( I ) *EQ.2) 82(1) «* 0 

COnTINUF- 

D0 300 I-1#160 

PRINT 4C0> I#R( I ) /Al ( I )#82( I ) 

FORMATnlO/lB/lB/Il) 

continue 

RETURN 

END 
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i: 

2 ! 

3: 

4 ; 
5J 
6 : 
7: 
8 : 
9 i 
10 ; 
11 ? 
12 ; 
13 : 
U; 
15: 
16 : 
17: 
18 ; 
19: 
20 : 
21 : 
22; 
23: 
24 ; 
25; 
26 : 
27; 
28: 
29; 
30; 
31 : 
32: 
33: 
34: 
35; 
36: 
37: 
38; 
39: 
40: 
4i; 
43 ; 
43: 


c 

C 

C 

C 


c 

c 

c 


600 


c 

G 

C 


44 

45 

46 

47 

48 

49 

50 

51 
5 ? 

53 


650 


R75-34 


VITERBI DECODING 

SUB'50UTr^‘E VDECaD(R) 

IMPLICIT InTEGER(A*Z) 

C0MM0N Ai (64O)/A?{64O)/PS{64O)iOR(64O)#R0(64O) 
dimension r(X) 

This pr0gram is for the viterbi DEceoiwG of a 
C0NV0LUTiaNAL C0DE 0F BINARY DIGITS/ WMERE K.3/ 

B»l/ V»2 ARE THE C6DE PARAMETERS 

dimension BM(4# 160)#ND(4/ 160 ) /RD ( 160 ) # DBM ( 4/ 160) 

INTEGER Bm 
integer DBm 
INTEGER HO 
integer R8 
integer RD 

N IS the desired length 9F THE BIT MEM0RY PATH PER PATH 

read ( 5#6 oO) N 
F0RMAT(I3) 

initialization 0 f decoding 

HD(1#1) * HD(2>D 9 HD(3/1 ) » mD( 4,1) » 0 

HD(t/2) « MD(2/2) • HDO/2) * hD<4^2) » 0 

SUMl » Al(l) + 1 

IF (SUMl .EQ,2) SVMl c 0 

SUM? « a2(D + 1 

IF(SUM2.EQ.2 ) sum? = 0 

SUM3 » Al{2) 4 1 

IF(SUM3 *Eq. 2) SUM3 - 0 

SUM4 * A2(2) + 1 

IF(SUM4.EQ,2) SUM4 - 0 

SUM5 • Al<3) 1 

IF(SUM5*CQ.2> SUM5 = 0 

SUM6 « a2(3) 4 1 

IF(SUM6*EQ,2) SUM6 a 0 

HDl a Al(l) + A2 (l> ♦ AK 2) + A2<2) + Al( 3) + A?(3) 

HO? a SUMl ♦ SUM2 + SUM3 + A?(?) SUMS 4 Sl'M6 
H03 a Al<l) + A2(l) + Al(2) 4 a2(2) + SUM5 4 SUM6 

HD4 a SUMl + SUM2 + SUM3 + A?(?) ♦ Ai(3) ♦ A?(3) 

HD5 • A1<1) + A2(l) + SUM3 SUH4 4 SUMS ♦ A2(3) 

HD6 a SUMl 4 SUM? 4 Al<2) 4 SUM4 + AK3) ♦ SUM6 

HD7 a Aid) + A?(l) ♦ SUMS 4 SUM4 4 AK3) SUM6 
HD8 a SIJMI 4 SUM? 4 Al(2) 4 SUM4 + SUM5 4 A?<3) 

PRINT 650 #hD 1/HD2/HD3#HD4/HD5/MD6/ND7/HD8 
F0RMAT(8l5) 

HD(1/3) » HD2 

BM(1 /I) thru BM( 4 /I) represent THE BiT MEM0RY PATHS 
AND HD(lil) thru H 0(4/D REPRESENT THf C0RRESP0NDING 

Hamming distance paths 

BM(l.l) ' 1 

ORIGINAL PAGE IS D-8 

OF POOR QUALITY 
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54 * 

55: 

56: 

57 ; 

5 tt j 

59: 

60 : 

61 ; 

6?: 

63: 

64 ; 

65': 

66 : 

67; 

68 : 

69 ': 

70 : 

7i: 

73 ; c 
73 ; 660 

74: c 
75: c 
76: c 
77- 
78: 

79: 

«0; 

?l : 

82: 

83: 

84 : 

85: 

86 ; 


wo i 

89 : 

9c; 

9i; c 

92; 670 

93: C 
94; C 
95; C 
96: 

97; 

98 : 

99: 

100 '- 

ici: 

1C2;' 

1C3; 

104: 


BM{1#2) ^ 0 

IF(HD1 •L^*HD2) HD(l#3)»HDliBM(lil) • 0j9M(l,2) « • 0 

HD (2/ 3) « HD4 

BM(2#n ^ 1 
BM(2#2) * 0 

BM{ 2 ^ 3 ) * 1 , * 

IF ( HD3#L£*HD4 ) HD(2#3)aH03i BM(2#l) »0| BM(?j2) 0;BM{2^3) • 1 

HD(3/3) “ HD6 
« 1 

BM(3i2) » I 
BM{3,3) » 0 

IF{MD5*LE*HD6) HD( 3#3) 3i 1) ■ 0;BM(3#2) • lJBM(3j3) > 0 
HD(4j 3) “ MD8 
BM(4,n • 1 
BM(4j2) ^ 1 

BM(4,3) = r , 

IF ( MD7#LE^«HD8 ) HD( 4# 3) eH07i bM( 4# 1) » 0;BH(4j2) ■ l;BM<4/3) * 1 
PRINT 660 #hD( l#3)#HD(2i 3) /HD(3/3 )>hD(4>3) 

F8RMAT(4l5) 

» * #'♦ » * ♦ ♦ ♦ ♦ • * 

D0 to I»4,160 

t ■ I • 1 

SUMl • Al(I) 4 1 

IF(<^UM1.F0*2) SUMl » 0 

SUN? « A2<I) 4 1 

IF{SUM2#EQ.2) SUM2 » 0 

HDl ¥ HD( 1/L) 4 Aid) 4 A?( I) 

HD2 a HDI3,L) 4 SUMl 4 SUMS 
HD3 a.HDdfL) 4 SUMl 4 SUMS 
HD4 = NDO-l) 4 Aid ) + A2d ) 

HD5 * HD<2iLi 4 SUMl 4 a2(I) 

HD6 « MD<4,t^ * Aid) 4 SUM2 
HD7 « HD(2#L) 4 A1(I) ♦ SUM2 
MD8 a MD(4^L> + SUMl 4 A2d) 

PRINT 670MiHDl/HD2iHD3/HD4/HD5/H06/HD7#HD8 
FeRMAT( 9l5) 

DBM<1#J) IS A dummy BIT MEMQRY PATH MATRIX 

BMdd ) -^ 0 
BM{2/I) « 1 
BM(3# 1 ) » 0 
BM<4#I) ^ 1 
Ml a 1 
M2 • I 

IF(I-6T*N) Ml § I • N + 1 
D6 20 MaMi,M2 
DBM(1>M) * BM(1>M) 
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1C5! 





106; 





107; 


■ BM(4^M) 



10«‘; 

20 

C0NTINUF 



109; 


M2 « I-*l 



no; 


IF(HD1.L5.h02) G» T0 672 


lit; 


HD(1jI ) * HD2 



11?; 


00 671 



113; 


BM(IjJ) » DBM(3iJ) 



lU: 

671 

C0NTINUF 



115; 


G0 T0 673 



116; 

672 

HDa#I) » HDl 



117; 

673 

IF(WD3*L5*MD4) G0 T0 675 


118; 


? MD4 



119! 


00 674' J»Mi,M2 



l2o; 


BM{2#J) « DBM(3#J) 



12l; 

674 

C0NTINUF 



1??; 


G0 T0 677 



123; 

675 

HD(2fl) ■ HD3 



124; 


D0 676 J»M1#M2 



125; 


BM(2#J) • DBM(UJ) 



126; 

676 

C0NTINUE 



127; 

677 

IF(HD5#U5«HD6) G0 T0 679 


128* 


H0(3#I) * H06 



129: 


Oe 678 J»M1#M2 



130'; 


BM(3#J) » DBM(4/J) 



131; 

678 

continue 



132: 


GO TR 682 



133; 

679 

HD(3# I ) ’ HD5 



134* 


06 681 J-M1/M2 



135: 


BM(3#J) « 0BM(2#J) 



136; 

681 

CONTINUE 



137: 

682 

IF(MD7fL5«H08) GQ T8 684 


138; 


HD(4W) » HD8 



139; 

683 

CONTINUE 



140: 


GO TO 685 



141 : 

684 

HD(4#n * HD7 



n?; 


DO 685 J^M1#M2 



1431 


BM(4i J) » DBM(2# J) 



144 ; 

685 

continue 



145: 


IFdfLT.N) GO TO 10 



n6; 

C 




147: 

C 

determines OLDEST 

BIT 0N 

THE 

148* 

C 




1^9; 


J 1 I . 



1-0 : 


K ■ I • N ♦ 1 



151 ; 


IF(HD( 1# J) ♦L2*HP(2/ J) ) 

GO TO 

740 

152'; 

760 

IF(H0(3/ J) •LE*HD(4# J) ) 

GO TO 

765 

153* 


IF(H0(2/ J) fLF»MD(4, J) ) 

GO TO 

710 

154 ; 


60 Te 730 



155; 

765 

IF(HD(2#J) .LE«HD(3/J) ) 

GO TO 

710 


M6ST path 
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156; Ge ra 7 ?o 

157; 740 IF (ND{ 3# J) .Lf •HD< J) ) GO T 6 770 

15?; IF(HDd^J) J) ) G9 TP 7C0 

159: G 6 7 ? 73 O 

160; 770 IF(HD{ 1 # J) J) ) G 6 T0 700 

I 61 ; G0 T 0 7?0 

16?; C 

163; C RP(T) IS decode 

164; C 

165; 700 R 8 (K) « BM( 1 #K ) 

166: G 8 re 800 

167; 710 R9{k') » BM(2>K) 

168; G 6 to 8 C 0 

169; 720 Re(k') * BM(3/K) 

170: 68 T^ 8 OO 

171; 730 R 8 (K) » Bm(4#K) 

17?; 800 continue 

173: 10 continue 
174: C 00 696 I’lil60 

175; r PRINT 695, I,BM( 1 ^ I )#3M(2j I )#P^(3# I )/BM(4# I 1 

176: r PRINT 695, I ,HD( Id ) #HD{ 2 # I ) #H0( I )#H0(4d ) 

177; 695 F8PNAT(5l5) 

178: 696 CBnTTNUE 

179 ! 06 950 J^ld^O 

lao; RD( I) » f^ 0 ( I ) ^ R( n 

181 ; PRINT 900 , I,R 8 ( I ),R( I),RD( I) 

1^2; 900 F0PMAT( Tl0,3l5) 

I83: 950 C 0 NTINUE 

184; RETURN 

185: end 


original pa§b is 

OF POOR QUALITY 
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LEVEL I ENCODE 


1 

? 

3 

if 

5 ; 

6 

7 

8 
9 

10 

11 

12 

13 

U 

15 


35 

36 

37 

38 

39 
^0 
^1 

4 if 


300 


4 C 0 


SUB^^QUTINE ENC0DE(NRITS/R ) 

IhPtlClT integer (A-Z) 
dimension R{1) 

CewMaN Al (6^0)/A2(6if0)#PS(640)/DRi6if0)#R0(640) 

DR(1 J«DR<3)»DR(3)»0 

oe 300 T»if#NBlTS’*’3 

DRn)«R<l- 3 ) 

continue 

DO 400 I«1/NBITS*3 
R( I )f0R( I ) 

continue 

PS(1)»1 

DO 1 I*1^N0ITS 

j« I + i 


K*I 42 


16: 

IE(R(I) 

• Eq 

.0 

• AND»R( J) 

•EQ-O.AND 

17: 

IF(R(I) 

• Eq 

.0 

•AND.R( J) 

•EQ.O. AND 

18; 

IE(R( I ) 

. EQ 

.0 

•AND*R( j) 

•EQ.O. AND 

19; 

IF{P(I) 

• Eq 

.0 

•AND*R( J) 

•EQ.O. AND 

20; 

IE(R( I) 

• Eq 

.0 

•AND.R(J) 

•EQ.O .AND 

2i; 

IF(R( I ) 

.EQ 

.0 

• AND.R(vl) 

. EQ.O. and 

22; 

IF(R(I) 

• Eq 

.0 

•AND*R( J) 

♦EQ.OfANO 

23: 

IE(Rn ) 

• Eq 

.0 

• AND«R( J) 

•EQ.O. AND 

24; 

IF(P(I) 

• Eq 

.0 

• AN0*R( J) 

•EQ.OfANO 

25: 

IF(R(I) 

• Eq 

.0 

• AnD«R( J) 

•EQ.O. AND 

26: 

IF(R( n 

• EQ 

.0 

•AND*R< J) 

•EQ. l.AND 

27: 

IF(R{ I ) 

• Eq 

.0 

•And ♦ r ( j ) 

•EQ. l.AND 

28; 

IF(R(I) 

.Eq 

.0 

•AND*R( J) 

•EQ«lf AND 

29; 

IF(R(I) 

.Eq 

.1 

♦AND*R{ J) 

•EQ^O.AND 

30: 

1F(R( I ) 

• Eq 

.1 

♦ AND • R { J ) 

•EQ-O.AND 

31 ■; 

IF(R( I ) 

•eq 

• 1 

•AND»R( J) 

•EQ.O. AND 

32: 

IF(R{ I ) 

« Eq 

.1 

tAND»R ( J) 

•EQ*1.AND 

33; 

IF<R( I) 

.Eq 

.1 

•AnD«R( J) 

•EO^l.ANO 

34'; 

IF{R(1) 

• Eq 

.1 

•AND*R( J) 

•EQ.l .AND 


AX( n«oj 

GO TO 1 
A1(I)-1J 
GO TO 1 
Al(t>*Cl 
GO TO 1 

AKn*o; 

GO TO 1 
Al(I)«0i 
GO TO 1 


A2(I)?0i PS(J)=1 
A2(l)»0i PS(J)-2 
A2(n«l^ PS(J)»3 
AH(I)»0; PS(J)«4 
A2(I)-l; PS(J>«4 


EQ.O 

EQtO 

EQ*0 


R(K) 

R(K) 

P(K) 

R(K) .EQ*0 
P{K) .EQ.O 
R(K) 

R(K) 

R(K) 

R(K) 

R(K) 

-PS( I ) .EO. 
PS( I ) .EO. 
.PS( I ) *E0* 
>PS(I).E0. 
.PS( I ) .EQ. 
,PS( I ) .EQ. 
■ PSd >.EQ. 
.PS( I ) .EQ. 
.PS( I ) .EQ. 


• AnD*PS( I ) 
.AnD*PS( I ) 
.AND«PS( I ) 
*AN 0 .PS( I ) 
.AND.PSt I ) 
EG.1.AND»PS( I ) 
EQ.1.AnD*PS( I) 
EQ.l .AND*PS( I ) 
EQ.1 .And«PS(I) 
EQ.l .A mD'PS ( I ) 


EQ.6) 
E0.3) 
EQ.l ) 
EQ.7) 
EQ.2) 
EQ.6) 
EQ.l) 
EO*2) 
EQ.3) 
E0.7) 


GO TO 
GO TO 
GO TO 
GO TO 
GO TO 
GO TO 
GO TO 
GO TO 
GO TO 
GO TO 


GO 

TO 

7 

GO 

TO 

7 

GO 

TO 

7 

GO 

TO 

8 

GO 

TO 

8 

GO 

TO 

9 

GO 

TO 

10 

GO 

TO 

10 

GO 

TO 

11 
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45; 

7 

A1(I)*»0; A 2 (I)pO^ pS(J>-5 

46; 


GO Ta 1 

47: 

8 

A1(I)»CJ A 2 (I)ai; PS(si )»6 

4» ; 


G 8 T 8 1 

49) 

9 

Al(T)»0J A2(I)»0I PS(J)»7 

50 : 


Ge T 8 1 

51: 

10 

AKDflJ A2(I)«0^ PS(J )«8 

5?: 


G 6 Ta 1 

53; 

11 

AKDpCJ A2(I)«0i PS(J)«9 

54; 

1 

caNTINUE: 

55: 


WRITE(6V172) 

56 ; 

172 

FaRMAT(lHl) 

57: 


D 8 80 I^1,NBITS/100+1 

58; 


IX«( I*l)*100+l 

59;, 


I2f Il-f99 

60: 


I2»miN{ I U99#NBITS) 

61 ; 


WRITE(6j85) ( R( J)/ J«I 1 # 12 ) 

6 ?: 


write ( 6 # 85) (AK J)^ J»11j I2) 

63 : 


WRITE (6^85) ( A2< J)/ J*I1> 12) 

64 : 

85 

F0RMAT(» SlOOIl) 

65 : 

80 

CeNTINUE 

66’; 

12 

CeNTINUF 

67; 


■ RETURN 

68 ; 


end 
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LEVEL I DECODE 


1 

2 

3 

4 

5 
6; 

7 

8 

9 

10 

u 

12 

13 

U 

15 ' 

U 

17 

18 
19 

ac' 

ai 

2 ? 

23 

24 

25 

26 

27 

28 

29 

30 

31 
3 ? 

33 

34 

35 

36 

37 

38 
39 : 
40 
^1 

42 

43 

44 


200 


subroutine 0 ECaDE(NBITS#R ) 

IMPLICIT INTEGER (A-Z) 

COMMON AK640) jA2(640)#PS(640)#DR(640)#R0(640) 
dimension bM( 9# 400)# HD (9/ 400)^ RO ( 400 ) # DB^ ( 9# 400 ) 

dimension B1<400)#B2(400) 
dimension rU ) 

Ni»l6 

EDIT«0 

dimension P( 200) #PP(200)#PM{9>200) 

P(D»P(2)«0 

F8RMAT{I3) 

HD(1# 1 ) »HD(2/1 ) «HD(3#1)»HD(4#1 )«H0(5/l )«HD(6# 1 ) *HD ( 7# 1 ) «HD ( 8# I ) -HO 
1 ( 9 / 1 ) «0 

HD(1/2) «HD{2/2) *HD ( 3# 2 ) »HD ( 4# 2 ) »HD (5# 2 ) *HD ( 6# 2 ) *^D ( 7# 2 ) »M0 ( 8# 2 ) pHO 

1(9#?)«C 

HD(?i3) «HD<3i3)«H0(7#3)«H0(9#3)»50 
SUM1»A1(1)+1 


IF(SIJM1,Eq, 2) 
SUM2*A2(1)+1 
IF(SUM2fEQ.2) 
SUM3*A1 (2)+l 
IF(SUM3.EQ,2) 
SUM4.A2(2)+1 
IF (SUM4.EQ.2) 
SUM 5 eAi( 3 )+l 
IF($UM5.E:q,?) 
SUM6»A2(3)+1 
lF(SUM 6 fE^Q* 2 ) 
HD(1#3) *SUM1 
HD(l# 3 )»SUMi 
HD(4#3) «SUM1 
HO { 5# 3 ) sSUmI 


SUMl.C 

SUM2«0 

SUM 3«0 

SUM4bO 


SUM5»0 


+A 1 ( 3 )+A 2 { 3 ) 
+ A 1 ( 3 )**'A 2 ( 3 ) 
+A 1 ( 3 )+A?( 3 ) 
+A 1 ( 3 )+A 2 ( 3 ) 


SUM 6*0 

■►A2( 1 )+A1(2)4.SUM4 
♦A2(l)+Al(2)+SgM4 
■►A2(1) + A1(2)+SUM4 

. _ +A2(1) + A1(2)^'SUM4 

H0(6#3) - Ak 1 )+A2( 1 )+A1 (2)+A2(2)+A1(3)-».SUM6 
HD(8#3)-Aia>4A2(l)+AX(2)+A2(2)+SUM5 +A2<3) 
P(3)*MIN0(HD< 1# 3)iHD(4#3)/HD<5#3)#Hn{6#3)#MD(8/3M 

BM{1#1)«9M(1#2)»BM(1/3)»0 
BM( 4# 1 ) sBm ( 4# 2 ) *BM ( 4#3 ) »0 
BM(5#l)«BM(5i2)»0; BM(5#3)«1 
BM(6# 1 ) ■BM(6/2) *0J BM{6 j 3)»1 
BM{«#l)-0; BM(8#2) sBM(8#3) ■! 
BM(?#l)«SM{?^2)eBM(2#3)=0 

BM ( 3 # 1 ) »Bm ( 3 # 2 ) »BM ( 3 # 3 ) .0 f 

BM( 7 , 1 ).BM( 7 # 2 )»BM( 7 # 3)«0 ^ 

BM(9#1)«BM(9#2)»BM(9j3)»0 ^ ^ 
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^ 6 : 
47 : 
48 5 
49 : 

so: 

5 ?; 
53 : 
54 : 
55 : 
56 : 
57 ': 
58 ; 
59 : 
60 : 
61 : 
6 ? ; 
63 : 
64 ; 
65 : 
66 ; 
67 : 


70 : 
71 ; 
73 ; 
73 : 
74 ; 
75 : 
76 ; 
77 ; 
78 ; 
79 : 
So: 
Si : 

8 a; 

83 : 
84 : 
85 : 
86 : 
87 : 
88 : 
89 : 
50 : 
3 i: 
92 : 
93 ; 
94 : 
95 : 


Oa 14 lp 4 ,NBlTS 

SUM 1 *A 1 (I )+1 
IF(SUM 1 ,EQ, 2 ) SUMl?iO 
SUM 2 »A 2 (I )+1 
IF(SLIM 2 «E:q* 3 ) SUM 2 iO 
HDl « H 0 ^ 6 /lp)+Al( I)+A 2 ( I) 

MD? « HD( 3 #C>+A 1 ( I )+A 2 ( I) 

HD 3 « HD( 1 j»L)+SUM 1 +A 2 ( I) 

HD 4 ■ ND< 7 ,L)^'SUM 1 +A 2 U) 

H 05 « HDt 2 #L>+Al(n+SUM 2 
HD 6 ■ HD< 6 jt^+Al ( I )+A 2 ( I ) 

HD 7 c WD<UL)+A 1 ( I )+A 2 ( I ) 

HD 8 ■ HD( 2 #t)’»'Al(I )+SUM 2 
HD 9 P yiO(3 ,l)*AX ( I )+A 2 ( I ) 

HP 10 » HD( 7 , 1 ^)+A 1 { ! )+A 2 ( I) 

HDll* HD< 4 it >+Al( I )+A 2 ( n 
HDl 2 « )-^Al ( I) + A 2 ( I) 

HDl 3 «P MD( 7 /L)+Al(I)+A 2 (n 
HD 14 - HDt 5 #t)+AlU )+SUM 2 
HDl 5 f WD< 9 #L)^A 1 U )«*'SUM 2 
HD 16 « Mn( 8 #L)+AK I)+A 2 ( I ) 

HD 17 « HD< 5 ,L>+SUM 1 +A 2 <I) 

HDl 8 » Mn( 9 ,l^)+SUM 1 +A 2 (I) 

HD 19 « HD< 8 #L)+AK I )+A 2 ( I) 

IF {EDIT*EQ*D 

•PRINT HDI^HD2^HD3^ HD4#H05i H06#HD7/H08#HD9#HDlOjHDll#H^12#HD13iH 

10 l 4 #M 0 l 5 iiHDl 6 /HD 17 ,HDl 8 jH 0 l 9 

15 F 8 RHAT( 20 I 4 ) 

QpMIN 0 (HD 1 ,HD 2 ) 

Ml.l 

M 2 *L 

IF(I.GT*N) mI.I-N +1 

D 0 16 MxMi,m 2 
DBm( 1 /M ) *BH{ 1 /M) 

DBM( 3 #M )«BM( 3 #M) 

DBM( 4 #M)*BH( 4 /M) 

DBm < 5 # M ) «BK ( 5 # M ) ORIGINAL PAGE IS 

DBm ( 6 # M ) » 8 m ( 6 # M ) OF POOR QUALITY 

DBM( 7 /M)*BM( 7 /M) 

DBM(8iM ) =BM(8#h) 

DBM(9/M)«BM(9/M) 

16 CdNTINUF 

IF (EDIT #NE* 1 ) G 6 T 0 l 7 l 
Oe 17 NlDXl«l /9 

PRINT 18 ^ I# (BM(NDX 1 #N 0 X 2 ) # NDXBxMl^Ma) 

18 F 0 RMAT( I 40 j 10 X# 19 I 3 ) 

17 C 8 NTINUE 
l 7 l C 0 NTINUE 
M2p I • 1 
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96 ; 

97 : 

SB: 

99; 

100 : 

ICl; 

102 ; 

103 ; 

104; 

105 : 

IC 6 : 

107 * 

1 C 8 ! 

109 '; 

110 : 

111 : 

112; 

113 ? 

114 ; 

115 : 

116 ; 

117 ; 

118 : 

119 : 

lao: 

131 ; 

132 : 

133 ’: 

134 : 

135 : 

126 : 

137 ; 

128 ; 

139 ; 

130 ; 

131 : 

132 ; 

133 : 

134 ; 

135 : 

136 ; 

137 ': 

138 : 

139 : 

140 : 

141 ; 

142 ; 

143 ; 

144 i 

145 : 

146 ; 


U*" I 

IF(Q«EGfWOl) G8 T8 19 
IF(Q#EQ#HD2) G8 Te 30 

19 HD(l#L)«»HDi; BMd^U^O 
00 38 

38 C0NTINUF 
60 T9 57 

20 NDdiD-^'Da; BMd/L)*0 
D9 39 Vf^dH2 
BMd#V)«0BM(3#V) 

39 CONTINUE 

57 QpHIN0(H03#HD4) 

IF(0.EQ*HD3) G0 T6 21 

if(q.eq.hd 4 ) g 0 re 22 

21 H0(2#U?HD3i BM(2#L)^0 
00 40 V»M1 jM2 
BM(2/V)»DBm(1#V) 

40. continue 
GO T9 58 

22 HD(2/L)«^^D4j BM(2^L)»0 
00 41 VpMi^M 2 
BM(2#V)»DBM(7iV) 

41 continue 

58 Q*MIN0(HD6#MD7#HD8#HD9/HD1C) 
HD(3#L)»HD5j BM(3#L)«»0 

DO 42 V«'^1#M2 
BM(3# V) *0bm(2# V) 

42 CONTINUE 
IF(0«EG*HD6) GB to 24 
IF(QtEG*H07) GO T0 25 
lF(0fEQ.H08) GO T6 26 
IF(Q.EQ*HD9) go to 27 
IF(Q*EQ«HD10) go to 28 

24 HD(4#U)*HD6; BM(4/L)^0 
DO 43 

BM{ 4 jV)» 0 Bl^( 6 #V) 

43 continue 
GO to 59 

25 HD(4^L)«H07; BM(4>L)»0 
Oe 44 V«Mi,n2 

BM(4# V) oOBMd/V) 

44 continue 
GO TO 59 

26 HD(4#L)«HD8J BM(4#L)«0 
DO 45 V«M1/M2 

8M(4# V) ■DBN<2/V> 

45 CONTINUE 
GO TO 59 

27 H0(4#L)pHD9I BM{4jL)«0 
DO 46 V*Mi^M2 
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147'i 


BM(4#V)«DBM(3#V) 

1<*8; 


CONTINUE 

l'*9; 


GO T8 59 

150; 

28 

H0(4jL)*HD 10^ BM{4iL)»0 

151: 


00 47 V«M1^M2 

152; 


BN(4#V)fOBM(7/V) 

153: 

47 

continue 

154; 

59 

Q»MIN0(HDU/H012#HD13) 

155'; 


IF(Q*EQ*HDn> GO T0 29 

156; 


IF(G).EQ.NDi2) G8 T0 30 

157; 


IF(0.EQ.HD13> go T0 31 

15g; 

29 

HD(5,L)-W0lli 8M(5#L)«l 

159'; 


DO 48 VfMl,M2 

160: 


BM(5iV)*GBM(4/V ) 

161: 

48 

continue 

162; 


GO TO 60 

163; 

30 

MD(5#U-HDl2i BM(5iL)»l 

164; 


DO 49 VtMl,M2 

165; 


BM(5>V)«Dbm(6/V) 

166': 

49 

continue 

167; 


GO TO 60 

168; 

31 

H0(5#L)«HD13I BN<5#L)«1 

169; 


DO 56 V«M1^N2 

170: 


Bli(5iV) »DBM(7/V) 

171; 

56 

continue 

172; 

60 

QfMlN0(HDl4/HDl5) 

173; 


IF<0*EQ«H014) go to 32 

174': 


IF(0*EQ*H015) GO TO 33 

175: 

32 

HD{6#U»Noi 4J BM{6iL)«0 

176: 


DO 50 V«Mv^M2 

177; 


8M(6/V)»DBH(5#V) 

178: 

50 

continue 

179; 


GO TO 61 

180: 

33 

HD(6#L)«HDl5i BM(6,L)-0 

I8i: 


DO 51 ViMl,M2 

182: 


BM(6#V) ■0bm(9# V) 

183; 

51 

continue 

184; 

61 

Q«NlN0(HDl7iHDl8) 

185; 


HD(7#L)fH0l,6l BM(7^L)»0 

186: 


DO 52 V«^liM2 

187'; 


BM(7# V) «DBm(8#V) 

188; 

52 

continue 

189; 


JF(Q*EQ*HD17) go to 35 

190; 


IF(QtE0*HDj8) GO TO 36 

191'; 

35 

HD(8#UfHDl74 BhM8jL)*l 

192; 


DO 53 VfMl,M2 

193; 


BM(8iV)fOBM(5iV) 

194'; 

53 

CONTINUE 

195; 


GO TO 37 

196; 

36 

HD{8#L)^HDl8i Bh(8iU)»l 

197; 


DO 54 V«M1,m2 
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I5g; 

199; 

2C0: 

aoi: 
202 ; 
2 C 3 ; 
204 ; 
205; 
206 : 
207: 
208; 
209; 
210 ; 
211 ; 
212 : 
213': 
2U; 
215; 
216 ; 
217 : 
218: 
219; 
220 ; 
221; 
222 ; 
223: 
224; 
225: 
226* 
227: 
228 ; 
229; 
230: 
23l! 
232; 
233: 
234; 
235; 
236'; 
237; 
238; 
239; 
240; 
241; 
242; 
243': 
244; 
245': 
246; 
247; 
248; 


BM(8#V)«0BM(9>V) 

54 CONTINUE 

37 HD(9#L)p^^0l9^ BM(9#U)»1 

DO 55 V»Mi,M2 
BM(9#V)»DBM(8#V) 

55 continue 
IF(T.LT»N) go to 14 
J»l 

K»I**N+1 

Q»rilN0(MD( J)#HD(2/ J)/HD(3> J) #HD{4iU)#H0(5# J)#H0(6> J)^H0(7# 

18 #J)#H 0 ( 9 # J) ) 

IF(0.E0.HD(1^J) ) GO TO 62 
IF{0.E0*WD^2/J) > Q8 TO 63 
IF(Q*EQ.HD(3#J) ) Q0 TO 64 
IF(QtEQ.WD{4i J) ) GO TO 65 
IF{QtEQ*H0(5# J) ) GO T9 66 
IF(Q.EQ*H0(6/J) ) GO TO 67 
IF{0*EQ.HD(7# J) ) G9 TO 68 
IF(Q*EQ*W0(8#J) > GO TO 69 
IF(QtEQ*HD(9# J) ) GO TO 70 

62 R0(K)»BM{l,K) 

GO TO 71 

63 R0{K)»BM(2#K) 

GO TO 71 

64 R0(K) f0M(3/K> 

GO T0 71 

65 R0<K)*BM(4#K) 

GO TO 71 

66 RO ( K ) »BH ( 5# K ) ORIGINAL PAGE IS 

GO TO 71 OF POOR QUALITV 

67 R0(K)«BM(6/K) 

GO TO 71 

68 R0(K)»BM<7#K) 

GO TO 71 

69 Re{K)«8M(8#K) 

GO TO 71 

70 R0(K)»B^^<9#K) 

71 continue 
14 continue 

IF (EDIT *NE* 1) 60 TO 721 

DO 72 leliNBiTS , ^ . 

PRINT 73^ I ,BM( 1# I )/BM(2/ I )iBM( 3i I) ^8M(4# I)#BN(5/ I I )iBM(7# I ) 

l#BN(8/ I )^BM(9/ n#HD( 1/ I )#HD(2# i )>HD(3> n#H0<4/ n#HD(5#I )#HD(6# I)#H 
2D ( 7# n # HD ( 8# I) i HD ( 9# 1) 

73 F0RMAT(19I5) 

72 CONTINUE 
721 continue 

WRlTe(6^l72) 

172 FOrHAT(IHI) 

00 76 IfliNBiTS 
R(I)iR(I+3> 
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349: 


Ren)»R8(I*2) 

350 : 

76 

continue 

35T: 


00 74 Iml,NBlTS 

352: 


RD(npR8(n*RU) 

353; 

74 

continue 

354: 


oe 741 1 «^ 1 #NBITS /100 + 1 

355: 


IU(I-1)*100*1 

356': 


I3?n>99 

357: 


I3fMIN(lU99^NBlTS) 

35, 8- 


WRjTE(6^85) (R0( J)# Jp 

359: 


WRlTE(6/^5) (R 

360: 


WRiTE(6#85) (RD< J)/U« 

36i: 

85 

F8RNAT(» SlOOIl) 

362: 

741 

continue 

363: 


RETURN 

364: 


END 
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